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infrared 0bse^tlons  have  indicated  the  strong  enhancement  of  certain  I 

created  in  ""h  ,  8  Pf ' ^  °f  aur°ral  activlty-  A  theoretical  model  has  been 

r  ^  !  dGf  StUdy  the  lnfrared  Processes  associated  with  the  aurora.  The  I 
311^^  3  i/e  r®sulJs  °f  this  study  depend  upon  a  thorough  knowledge  of  the  rates  of 

we 1  have  Lde JT  h  h  “  °f  ^  ^  C°nStantS  are  Wel1  kno-  while  *>r  other! 
indirar^h  ,  H  crudest  approximation.  Qualitatively,  however,  our  results 

during' Der ind  "T""  "  ^  bands  °f  C°^  ’  N0>  and  N0+  are  stiongly  enhanced 

during  periods  of  auroral  activity.  This  is  primarily  the  result  of  two  basic  energy 

xchange  me.  hamsms.  The  energy  deposited  by  the  primary  auroral  electron  flux  goes  ' 

»~rIS“’7„th:.ln,“rjS  enef8V  0£  N”  «*■  i-  the  vibrational  ! 

This  energy  is  then 5  significantly  increased  in  an  aurorally  disturbed  atmosphere.  j 

as  So!  I  X  NO  !tc  rh  ^  ^  Varl°US  infrared  ac^e  "inor  constituents  such  ! 

cesses’  increasing  ’  Jbr°ugh  nea[  reeonant  vibrational-vibrational  exchange  pro- 
tblr  V  Increasing  the  higher  level  vibrational  opulations  of  these  species  and 

pamc  e  t:;L Le‘r  Inf”red  r:diatiVi!  0utput-  In  ^“ion  to  energetic3 auroral  1 
the  results  are  simff  8°UrCe’  Wa  have  also  considered  electric  fields.  Qualitatively 
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NOTE  ADDED  IN  PROOF 


The  radiances  computed  for  selected  bands  of  the  infrared  active 
minor  constituents  treated  in  this  study  are  affected  by  various  input 
radiative  fluxes,  as  discussed  in  the  report.  The  values  adopted  for 
the  terrestrial  flux  were  based  on  the  study  of  Corbin  et  al.,  [1970]. 

It  has  recently  been  suggested  [ Murphy ,  1973],  and  subsequently 
confirmed  [Degges,  1974],  that  the  radiative  flux  from  the  lower 
boundary  for  the  present  cal.  dations  (60  km),  as  provided  by  the 
cited  study,  is  probably  significantly  in  error  over  portions  of  the 
spectral  region  of  interest.  In  strongly  absorbing  spectral  legions, 
the  effective  blackbody  temperature  at  the  lower  boundary  is  expected 
to  closely  approach  the  local  kinetic  temperature.  (For  example,  in 
the  spectral  region  around  4.3  microns,  where  the  atmosphere  is  op¬ 
tically  very  thick  at  60  km,  the  appropriate  blackbody  temperature 
should  approximate  250‘L  for  the  temperature  profile  assumed  in  the 

present  study,  whereas  a  value  of  220°K  was  obtained  from  the  work  of 
Corbin  et  al.) 

The  general  impact  of  the  error  is  to  render  too  small  the  cal¬ 
culated  vibrational  populations  (and  hence  band  radiances),  particu¬ 
larly  in  the  region  of  the  lower  boundary.  Corrections  must  be  made 
on  a  band-by-band  basis  since  the  local  optical  thickness  differs  from 
band  to  band.  The  significance  of  the  corrections  will  also  depend  on 
the  relative  importance  of  the  various  mechanisms  which  determine  the 
population  of  a  particular  excited  level. 

It  is  important  to  note  that  the  problem  with  the  terrestrial  flux 
does  not  invalidate  the  primary  conclusion  of  the  study,  i.e.,  that  the 
auroral  energy  input  to  the  upper  atmosphere  can  significantly  enhance 
certain  atmospheric  infrared  emissions  during  and  immediately  after 
periods  of  strong  auroral  activity. 

Among  all  the  bands  considered  in  this  study,  probably  the  CH4  band 
at  3.31  p  is  the  one  which  has  the  biggest  error  because  the  absorption 
of  terrestrial  flux  plays  the  predominant  role  in  excitation  of  this 
band.  An  approximate  estimation  indicates  that  differences  as  large 
as  one  order  of  magnitude  could  be  caused  by  the  error  in  terrestrial 
flux  on  the  calculated  band  radiance. 
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INTRODUCTION 


rhe  purpose  of  *iiis  study  is  to  investigate  the  natural  infrared 
background  or  the  earth  in  the  1  to  25  micrometer  region  of  the  electro 
magnetic  spectrum.  The  study  is  primarily  concerned  with  the  natural 
infrared  backgn  ind  during  aurorally  disturbed  conditions.  in  order  to 
define  this  condition,  the  quiescent  background  must  also  be  defined. 

It  was  originally  felt,  and  has  been  substantiated  by  this  study,  that 
the  auroral  energy  input  would  significantly  change  some  of  the  atmo¬ 
spheric  infrared  emissions  during  times  of  auroral  activity. 

The  emphasis  of  the  work  effort  was  the  development  of  a  computer 
program  to  physically  model  the  atmosphere  in  an  effort  to  predict 
infrared  emissions  under  various  conditions  of  auroral  activity.  As 
a  starting  point,  the  atmospheric  radiance  models  developed  by  Corbin 
et  2l.  [1970]  under  a  Honeywell,  Inc.  contract  let  by  the  Air  Force 
Cambridge  Research  Laboratories  were  adopted.  The  so-called  low-alti¬ 
tude  program  was  used  only  to  determine  terrestrial  fluxes  as  an  input 
to  the  new  program  to  be  developed.  The  so-called  high-altitude  pro¬ 
gram,  on  the  other  hand,  was  used  as  a  starting  point  for  the  new  pro¬ 
gram. 

Tae  high-altitude  programs  of  Corbin  et  at.  [1970]  predict  infra¬ 
red  emissions  for  the  gases  CH„,  C02 ,  H20,  NO,  N20,  and  03  during  qui¬ 
escent  conditions.  These  programs  contain  no  atmospheric  chemistry. 

The  present  study  undertook  to  modify  these  programs  by  .iding  addi¬ 
tional  minor  constituent  gases,  extending  the  spectral  region  to  in¬ 
clude  some  emissions  from  the  region  1.0  to  5.0  micrometers,  and  lastly 
but  more  importantly,  to  include  a  time-dependent  accounting  of  the 
effect  of  the  auroral  energy  input  on  the  vibrational  populations  as 
well  as  the  number  density  profiles  of  the  various  minor  constituents 
under  consideration.  This,  of  course,  implies  that  a  number  of  atmo¬ 
spheric  chemical  processes  were  considered. 

"'he  high-altitude  program  of  this  study  includes  infrared  emis- 
sions^from  CH4 ,  C02 ,  H20,  NO,  N20,  03,  N0+,  OH,  and  CO.  The  inclusion 
of  NO  in  the  list  of  gases  changes  the  context  of  minor  constituents 
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to  that  of  ions  and  neutral  molecules  rather  than  solely  neutral  mole¬ 
cules.  The  programs  which  account  for  atmospheric  chemistry,  therefore, 
must  not  only  be  concerned  with  neutral  chemistry  but  must  also  be  con¬ 
cerned  with  ion  chemistry,  as  well  as  the  interaction  of  charged  auroral 
particles  with  the  atmosphere.  A  major  programming  effort  was  also 
necessary  in  order  to  determine  the  effects  that  the  incoming  auroral 
particles  would  have  on  the  vibrational  temperature  of  molecular  nitro¬ 
gen.  An  accurate  knowledge  of  the  vibrational  population  of  molecular 
nitrogen  is  necessary  because  of  the  importance  of  the  many  near-reso¬ 
nant  vibrational -vibrational  energy  exchange  processes  between  N2  and 
minor  constituent  gases. 

One  of  the  results  of  the  current  study  is  in  the  form  of  a  plot 
or  a  table  giving  the  infrared  band  radiance  in  watts/cm2-steradian  for 
various  important  bands  of  the  minor  constituent  gases  previously  men¬ 
tioned.  The  infrared  radiance  is  given  in  5-km  increments  from  60  km  up 
to  15C  Km  (or  the  point  where  the  number  densities  of  the  gas  become  in¬ 
significant)  and  then  in  25-km  increments  thereafter  to  500  km.  The 
direction  in  which  the  observer  is  looking  is  varied  and  includes  the 
vertical  direction  (zenith  angle  0°)  as  well  as  zenith  angles  of  20°, 

40  ,  60°,  70°,  80°,  and  90°  (horizontal).  One  other  geometry  (limb 
view)  is  considered  and  that  is  me  placing  the  observer  outside  of 
the  atmosphere  and  allowing  him  to  view  along  a  path  which  is  tangent, 
not  to  the  earth’s  surface,  but  rather  to  an  atmospheric  shell  above 
the  earth.  The  atmospheric  shells  (points  of  tangency)  are  located  in 
5-km  increments  from  60  km  to  150  km  and  thereafter  in  25-km  increments. 
The  height  from  the  surface  of  the  earth  to  the  atmospheric  shell, 
which  is  tangent  to  the  view  path,  is  called  the  tangent  height.  The 
viewpath  of  the  observer  does  not  stop  here  but  rather  extends  to  a 
point  outside  of  the  atmosphere  on  the  side  away  from  the  observer. 

A  second  form  of  result  is  presented  as  a  plot  or  a  table  giving 
a  compilation  of  the  infrared  radiance  of  all  the  bands  of  all  the 
gases.  This  is  done  for  the  limb  view  case,  the  vertically  looking 
case,  and  the  horizontally  looking  case.  The  table  or  plot  lists  the 

radiance  in  watts/cm2-steradian  every  0.1  micrometer  over  the  region 
of  interest. 
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Several  auroral  and  quiescent  sets  of  conditions  were  studied. 

In  this  context  a  quiescent  atmosphere  has  the  meaning  that  no  vis¬ 
ible  aurora  is  present.  We  have  examined  the  infrared  output  in  the 
I-/5  micron  spectral  region  under  this  definition  of  quiescence  but 
with  an  electron  drizzle  present.  We  have  also  taken  a  quiescent  case 
but  with  electric  fields  being  present  in  the  atmosphere.  Further, 
we  have  studied  the  infrared  emissions  from  the  atmosphere  with  a 
class  IBC-II  and  a  class  IBC-III  aurora  in  progress.  A  more  detailed 
description  of  the  input  parameters  during  both  quiescent  and  aurorally 
disturbed  conditions  will  be  presented  in  a  following  section  of  this 
report . 


GENERAL  DESCRIPTION  OF  IR  BAND  MODEL  PROGRAM 

The  purpose  of  this  section  of  the  report  is  to  give  a  brief  but 
adequate  description  of  the  general  philosophy  used  in  computing  the 
Infrared  emissions  on  t.  band-by-band  basis  for  each  of  the  minor  con¬ 
stituent  gases  under  consideration. 

The  quantity  of  infrared  radiation  emitted  by  any  band  is  directly 
proportional  to  the  number  density  of  molecules  in  the  excited  state 
capable  of  radiating  at  the  band  wavelength  .  It  is  therefore  necessary 
to  have  as  an  input  to  the  model  the  best  information  available  with 
regard  to  the  height  profile  of  the  number  density  for  each  species 
under  consideration.  With  a  knowledge  of  the  height  profile  and  a 
given  temperature  model,  one  may  easily  calculate,  using  the  Boltzman 
equation,  the  populations  of  the  various  vibrational  levels  of  a  gas. 
These  results  represent  the  local  thermodynamic  equilibrium  (or  LTE) 
case.  At  auroral  altitudes  one  is  often  not  interested  in  this  case, 
since  many  atmospheric  conditions  do  not  justify  the  assumption  of 
lo.  al  thermodvnamic  equilibrium.  When  local  thermodynamic  equilibrium 
londiiions  do  not  prevail  (the  NLTE  case),  one  must  consider  radiative 
energy  exchange  processes  as  well  as  collisional  energy  exchange  pro¬ 
cesses  in  order  to  determine  with  reasonable  certainty  the  vibrational 
populations  of  a  gas. 
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In  order  to  compute  the  vibrational  populations  of  a  gas  under 
NLTE  conditions,  one  must  have  as  an  input  to  the  program  such  quanti¬ 
ties  as  solar  flux,  terrestrial  flux,  and  radiative  flux  from  other 
parts  of  the  atmosphere,  as  well  as  excitation  and  de-excitation  rates 
for  the  collision  processes  of  vibrational-translational  (V-T) ,  trans¬ 
lational-vibrational  (T-V) ,  and  near  resonant  vibrational-vibrational 
(V-V)  energ/  exchange.  Additional  inputs  giving  information  as  to  the 
characteristics  of  the  molecule  itself  are  also  required.  These  would 
include  the  energy  levels  and  their  degeneracies,  rotational  constants, 
band  transitions,  and  band  strengths. 

With  this  background  we  may  now  outline  the  methodology  used  in 
the  computation  of  the  emission  of  a  particular  infrared  band.  First 
the  required  inputs  as  described  above  are  read  into  the  computer.  A 
Boltzman  or  LTE  population  is  determined  for  each  of  the  vibrational 
levels.  Then,  using  this  population  as  a  starting  point,  an  equation 
describing  both  radiative  and  collisional  excitation  and  de-excitation 
processes  is  solved  in  an  iterative  fashion,  leading  to  a  new  value  of 
vibrational  population.  The  new  population  is  new  representative  of 
NLTE  conditions. 

Although  a  study  of  how  the  vibrational  populations  of  various  gas¬ 
es  change  with  changing  NLTE  conditions  is  of  interest  in  itself,  it  is 
more  important  to  relate  this  to  a  quantity  which  can  be  measuied.  We 
desire  to  know  the  band  radiance  expressed  in  this  report  as  watts/cm'- 
steradian  that  one  would  observe  for  various  infrared  bands  for  a  variety 
of  conditions.  If  one  integrates  this  quantity  over  a  column  with  unit 
area  cross  section  along  some  optical  path  length,  one  obtains  the  pre¬ 
dicted  radiance  for  the  chosen  band  and  NLTE  conditions.  Most  of  the 
gases  considered  are  optically  thin;  however,  in  the  case  of  the  4.3  p 
band  of  CO j  this  is  not  true.  In  this  case  self  absorption  must  be  in¬ 
cluded  in  the  integration  in  order  to  obtain  the  radiance. 

The  picture  which  has  been  presented  thus  far  appears  simple 
enough  on  the  surface;  however,  in  a  detailed  analysis,  it  becomes 
rather  complicated,  especially  during  periods  of  electron  precipita¬ 
tion.  Let  us  briefly  consider  the  procedures  involved  in  handling 
the  energy  exchange  processes  during  an  auroral  distrubance.  First  a 
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model  is  chosen  to  represent  the  primary  auroral  electron  flux.  With 
a  knowledge  of  the  primary  flux,  one  may  calculate  the  energy  deposi¬ 
tion  rate  for  the  electrons.  Bv  applying  the  proper  cross  sections, 
the  total  production  rates  of  Nj,  and  0"*"  may  also  be  calculated 
These  production  rates,  in  addition  to  various  neutral-neutral  and 
ion-neutral  chemical  reactions,  charge  exchange,  and  recombination 
must  be  considered  in  order  to  determine  the  height  profiles  for  such 
quantities  as  [n2],  [02],  [o  ],  [NO  ],  and  [NO].  These  height  pro¬ 
files  are  computed  on  a  time-dependent  basis  in  order  to  study  changes 
which  take  place  shortly  after  an  aurora  commences.  It  has  been  found 
that  of  the  minor  constituents  under  consideration  the  height  profiles 
of  [NO]  and  [N0+]  change  in  a  significant  enough  manner  after  auroral 
commencement  to  warrant  usin^  these  new  height  profiles  as  an  updated 
input  to  the  program.  In  addition  to  the  height  profiles  just  de¬ 
scribed,  a  production  rate  for  secondary  electrons  may  also  be  com¬ 
puted.  This  has  been  done  following  the  techniques  of  pees ,  Stewart, 
and  l Valuer  [1969].  With  an  accurate  knowledge  of  secondary  electron 
production,  one  may  compute  the  secondary  electron  flux  and  then  in 
turn  (and  with  a  knowledge  of  the  cross  sections)  the  production  rates 
for  N2(v=1),  02(v=1),  and  O^D). 

From  the  description  just  given,  it  can  be  seen  that  a  large  por¬ 
tion  of  the  auroral  energy  is  at  least  temporarily  deposited  in  the 
form  of  vibrational  energy  for  the  molecular  N2  and  02  and  in  the 
form  of  electronic  energy  for  atomic  oxygen.  At  this  point  several 
other  energy  exchange  processes  occur.  Energy  stored  in  the  0(JD) 
atom  becomes  an  important  production  source  for  N2(v=l)  while  both 
N2(v-1)  and  02(v=l)  exchange  energy  through  near  resonant  (V-V)  pro¬ 
cesses  with  all  of  the  minor  constituents  except  OH. 

Study  has  shown  that  due  to  the  importance  of  the  near  resonant 
(V-V)  energy  exchange  process ,  several  bands  are  particularly  ennanced 
during  auroral  activity.  These  include  the  bands  of  C02 ,  NO,  N20,  and 
N0+. 

Thus  it  is  seen  that  in  order  to  accurately  predict  infrared 
emissions  (particularly  for  the  bands  mentioned  above)  one  must  have 
available  an  accurate  knowledge  of  the  vibrational  population  of  molec- 
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ular  nitrogen.  The  methods  used  in  determining  the  vibrational  popu¬ 
lation  of  N2  will  be  presented  in  detail  in  a  later  section. 

INPUT  PARAMETERS  FOR  THE  INFRARED  BAND  MODELING  PROGRAM 

The  height  profiles  of  the  number  densiiies  of  the  various  minor 
constituent  gases  as  well  as  the  height  profiles  for  [o],  [o2l,  and 
[N2]  and  the  temperature  were  obtained  from  a  variety  of  sources.  This 
also  ho. -is  true  for  the  molecular  constants,  the  band  strengths,  the 
input  fluxes,  etc.  It  is  therefore  the  purpose  of  this  section  of  the 
report  to  collect  together  what  are  considered  the  more  important  of 
these  parameters  and  to  state  the  source  of  the  values  used. 

Number  Density  Height  Profiles 

The  height  profiles  of  [CII4],  [C02],  [h2o],  [n20],  [o3]  noon,  [03] 
predawn,  [oh]  day,  lOH]  night,  and  [CO]  are  shown  in  Figures  1-9  re¬ 
spectively.  The  height  profiles  of  [CH4],  [C02],  [h20],  and  [n20]  were 
adopted  from  sorbin  et  at.  [1970].  It  was  necessary  to  modify  the  [03] 
predawn  height  profiles  of  Corbin  et  at.  in  order  to  make  them  compat¬ 
ible  with  the  [OH]  profiles  of  Shimazaki  and  Laird  [1970].  To  explain 
this  further,  it  was  assumed  in  this  report  that  OH  is  predominantly 
formed  by  the  process  H  +  03  -*  OH  +  02 .  We  have  used  the  height  pro¬ 
files  of  both  [h]  and  [OH]  (day  and  night)  from  Shimazaki  and  Laird 

[l970].  The  [o3j  profile  used  by  Corbin  et  dl.  was  somewhat  larger  at 
higher  altitudes  than  that  of  Shimazaki  and  Laird  and  had  we  used  it, 

it  would  have  led  to  a  population  inversion  >f  [OHj  as  it  is  modeled 

in  this  program.  On  the  other  hand,  the  profiles  by  Shimazaki  and 
Laird  do  not  have  quite  the  structure  at  1>  ar  altitudes  as  those  of 
L^rbrn  er  at.  Therefore,  it  was  decided  to  simply  modify  the  height 
profiles  of  Corbin  et  at.  in  the  higher  altitude  region.  The  one  re¬ 
maining  height  profile  is  [CO].  It  was  assumed  that  CO  constitutes 
7 .0  x  10-®  times  the  total  number  density  as  suggested  by  Goody  [1964]. 

Figures  10-13  show  the  height  profiles  of  [NO]  and  [N0+]  during 
an  electron  drizzle  but  no  aurora  and  during  an  IBC-III  aurora  respec- 


Figure  1  .  Height  profile  of  CII4 


Figure  4.  Height  profile  of  N2O. 


mm 


profile  of  NO  ,  no  aurora  but  with  an  electron  drizzle  presert. 
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tively.  These  height  profiles  are  computed  in  the  chemistry  subprogram 
which  will  he  explained  in  detail  in  a  later  section  of  this  renort. 

The  height  profiles  of  [n2],  [O;],  and  [o]  are  shown  in  Figures  14, 
15,  and  16.  These  profiles  have  been  extracted  from  Adams  [1970].  The 
height  profile  of  the  total  number  density  is  shown  in  Figure  17. 

Figure  18  shows  the  temperature  model  used.  This  model  is  the  spring/ 
till  mid-latitude  model  from  the  U.S.  Standard  Atmosphere  Supplement 
[ 1966] . 


Rotational  Constants 

Of  the  molecules  under  consideration  in  this  report,  all  are  linear 
with  the  exception  of  CH^ ,  H20,  and  O3.  The  rotational  constants  Bv  of 
a  diatomic  molecule  for  a  given  vibrational  level  v  are  different  from 
the  rotational  constants  for  the  equilibrium  positions,  B  .  To  find  ^ 
an  averaging  process  is  used  which  can  be  shown  to  yield 

Bv  =  Be  -  ae(v  +  h)  +...  (1) 

where  oie  is  a  constant  which  depends  upon  the  anharmonicity  of  the 
vibrations  as  well  as  Be  and  other  factors.  For  linear  polyatomic 
molecules  the  expression  for  By  becomes 

r  d-i 

Bv  "  Be  -  la±  (vt  +  +.  .  .  (2) 

Here  the  subscript  v  has  become  the  set  of  quantum  numbers  vj ,  v2 , 

V3,...  and  dj  =  1  or  2  for  nondegenerate  and  doubly  degenerate  vibra¬ 
tions  respectively.  Table  1  is  a  compilation  of  the  rotational  con¬ 
stant,  Bv  for  all  of  the  minor  constituents  along  with  the  source  of 
these  constants.  Also  listed  in  this  table  are  the  wave  numbers  of 
the  vibrational  levels.  These  were  taken  either  from  the  listed  source 
or  CorHn  et  a' .  [1970] . 


Band  Strengths 

The  integrated  absorption  S^u  with  dimensions  cm-1  sec-1  may  be 
written  in  the  following  manner: 


NUMBER  DENSITY  (cm~3) 


NUMBER  DENSITY  (cm*3) 


NUMBER  DENSITY 


Figure  18.  Temperature  height  profile. 
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TABLE  1 

MOLECULAR  CONSTANTS 


Molecule 


CHU 


C02 


NO 

v  =  0 
v  =  1 
v  =  2 


vi  v2  v  3 


0  0  0  0 
0  10  0 
0  0  10 
0  0  0  1 


0  0  0 
10  0 
0  0  1 
0  l1  o 
0  2°  0 
0  22  0 
0  31  0 

0  33  0 


RCM 

Bv 

Source 

0000.0 

0.3400 

Corbin  et  al •• 

1526.0 

0.3500 

[1970] 

3018.4 

0.3500 

1306.2 

0.3500 

0000.0 

0.38918 

Hersberg  [ 1945] 

1388.2 

0.38862 

2348.0 

0.38628 

667.  < 

0.38180 

1285.4 

0.38990 

1335.1 

0.39042 

1932.5 

0.30104 

2003.3 

0.39104 

0000.0 

1.6961 

Albritton 

1878.0 

1.6785 

[1973] 

3727.0 

1.6610 

N20 


0  0  0 
l  0  0 

0  0  1 
0  1  o 
0  2  0 
0  3  0 

0  4  0 


0000.0  0.41892 

1285.0  0.41549 

2223.5  0.41984 

588.8  0.41916 

1173.0  0.41941 

1758.8  0.41966 

2338.5  0.41991 


Thomson  and  Wil¬ 
liams  [1953] 


0000.0  1.9887  Field  [1973] 

2344.0  1.9696 

4660.0  1.9508 


OH 

v  =  0 
v  =  1 
V  =  2 
v  =  3 
v  =  4 
v  =  5 
v  =  6 
v  =  7 
v  =  8 
v  =  9 


0000.0  18.513 

3562.0  17.805 

6976.0  17.097 

10222.0  16.389 

13315.0  15.681 

16260.0  14.973 

19061.0  14.265 

21725.0  13.557 

24258.0  12.849 

26663.0  12.141 


Herzberg  [1971] 
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TABLE  1  (cont.) 


Molecule 

CM 

n— ♦ 

> 

v3 

v4 

RCM 

B 

V 

Source 

CO 

v  =  0 

0000.0 

1.9255 

Mantz  et  al. 

v  =  1 

2143.0 

1.9050 

[1971] 

v  =  2 

4260.0 

1.8875 

h2o 

0 

0 

0 

0000.0 

0.3000 

Corbin  et  al. 

0 

1 

0 

1594.7 

0.3000 

[1970] 

03 

0 

0 

0 

0000.0 

0.3000 

Corbin  et  al. 

0 

0 

1 

1043.0 

0.3000 

[1970] 

Si  " 

lu 


Av 


lu 


K  d  =  (N  B  -  N  B  .) 
L,v  v  1  l-*-u  u  u+l' 


hv 


lu 


(3) 


where  k^v  (in  cm  is  the  spectral  absorption  coefficient  per  unit 

length,  and  where  (N^)  Is  the  number  density  in  the  upper  (lower) 

quantum  state  and  B^^  (B^^)  tlie  Einstein  coefficient  for  induced 
absorption  (emission).  The  parameter  is  the  frequency  (sec-1)  of 
the  band  center.  Often  it  is  advantagaov**  to  define  a  spectral  ab¬ 
sorption  coefficient  (in  cm-1,  atm'1)  where  uj  is  the  wave  number 
(u)  =  v/c)  .  It  is  the  integral  of  P  that  leads  to  the  band  strength 
Slu  as 


S 


lu 


P  dw 
u 


(4) 


In  this  expression,  p  is  the  partial  pressure  of  the  absorber  (in  atm) . 
Needless  to  say,  P^  (k^  v)  is  nonzero  only  for  a  narrow  wave  number 
(frequency)  range  (Avlu) . 

The  band  strengths  for  all  but  N0+,  OH,  and  CO  were  taken  from 
Corbin  et  al.  The  band  strengths  of  NO-1"  were  obtained  from  Stair  and 
jauvin  [1967].  Those  for  OH  and  CO  were  taken  from  Penner  [1959]. 

Table  2  lists  the  various  hands  of  interest  and  their  respective  band 
strengths . 
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fABLE  2 


9 

4.26 

H20 

1 

6.27 

5  x  IO-11 

2 

x  IO'15 

NO 

1 

5.32(1-0) 

4  x  10"15 

2 

x  10"15 

2 

2.68(2-0) 

4  x  10-15 

2 

x  IO"15 

3 

5.41(2-1) 

4  x  10"15 

2 

x  10-15 

n2o 

1 

16.98 

10-15 

2.5 

x  IO-I3 

2 

7.78 

10-15 

2.5 

x  10-13 

3 

4.50 

io-15 

2.5 

x  IO"!3 

°3 

1 

9.59 

10"14 

IO"34 

NO+ 

1 

4.27 

8  x  10-!5 

6 

X  10"15 

2 

2.15 

8  x  10"15 

6 

x  10"15 

3 

4.32 

1 

O 
•— < 

X 

00 

6 

x  10-i5 

OH 

1(7-5) 

1.87 

10"17 

io"i6 

2(4-2) 

1.58 

10-17 

10-16 

CO 

1 

4.67 

10"16 

10-15 

2 

2.35 

10"16 

10-15 

2700.0000  299 

300.0000  273 

125.0000  273 

2.1100  273 

0.0310  273 

33.0000  303 

260.0000  303 

1856.0000  303 

360.0000  273 

500.0000  300 

3.0000  300 

0.0100  300 

4.0000  300 

4.0000  300 
260,0000  300 
1.7900  300 
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Input  Radiative  Fluxes 

Radiative  input  fluxes  come  from  three  sources:  the  sun,  the  earth, 
and  other  parts  of  the  atmosphere.  In  the  computational  model,  the 
solar  and  the  terrestrial  fluxes  have  fixed  input  values  for  each  band, 
whereas  the  flux  from  other  parts  of  the  atmosphere  is  computed  in  the 
program  itself .  The  input  values  for  the  solar  flux  for  all  bands  ex¬ 
cept  NO  ,  OH,  and  CO  were  obtained  from  Corbin  et  at.  [1970],  whereas 
the  solar  flux  for  the  bands  of  these  gases  were  adapted  from  Gast 
[1965].  The  input  values  for  the  terrestrial  flux  were  also  obtained 
from  orbin  et  at.  for  all  but  the  same  three  molecules.  For  the  vi¬ 
bration-rotational  bands  of  these  molecules  the  so-called  "low  altitude 
computer  program"  [ Corbin  et  at.  1970]  was  run  and  values  obtained 
from  this  source.  For  band-center  wavelengths  below  the  operational 
limit  of  this  program  (5p)  ,  an  extrapolation  was  made.* 

Collisional  Energy  Exchange  Coefficients 

Four  basic  collisional  energy  exchange  processes  have  been  con¬ 
sidered  in  this  study.  They  are  vibrational-translational  (V-T) , 
translational-vibrational  (T-V) ,  vibrational-vibrational  (V-V) ,  and 
internal  vibrational  mode  exchange.  Of  these,  the  latter  Is  used  only 
for  C02  and  the  temperature  dependent  values  for  this  coefficient  come 
from  Taylor  and  Bitterman  [1969]. 

The  V-T  and  T-V  coefficients  become  important  to  the  excitation 
and  de-excitation  of  N2  and  02  (between  v  =  0  and  v  =  1  levels)  in  the 
upper  E-region  and  the  F-region.  The  V-T  coefficients  were  chosen  to 
be  an  input  parameter  to  the  program  and  the  T-V  coefficients  were 
found  by  solving  a  detailed  balance  problem  as  part  of  the  overall  com¬ 
puter  program.  The  V-f  values  used  by  Corbin  et  at.  were  used  for  all 
gases  except  NO  ,  CO,  and  )H.  No  information  could  be  found  regarding 
this  parameter  for  N0+.  As  a  result,  an  estimated  value  of  twice  the 
V-T  coefficient  for  NO  was  used.  The  same  problem  of  a  lack  of  avail¬ 
able  information  was  encountered  for  CO  and  OH.  This  resulted  in  es¬ 
timates  being  used  for  these  gases  also. 

The  reader  is  referred  to  the  Note  Added  in  Proof,  p.  x. 
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Of  the  collisional  energy  exchange  processes,  the  V-V  exchange 
process  is  probably  the  most  important.  This  is  because  near  resonant 
exchange  conditions  exist  between  N2  and  the  following  gases:  C02 ,  NO, 
N2C,  NO  ,  and  CO.  Also,  but  less  importantly,  near  resonant  exchange 
conditions  exist  between  02  and  CH„  ,  H20,  and  03.  Again  the  values 
used  by  ncrbin  et  al.  were  adopted  for  all  but  N0+  and  CO  (OH  having 
no  near  resonant  V-V  exchange  with  either  N2  or  02).  A  lack  of  informa¬ 
tion  was  again  a  problem  for  the  determination  of  an  input  value.  An 
estimated  value  of  three  times  the  NO  value  was  used  for  N0+  and  an 
estimate  consistent  with  the  other  constituents  was  used  for  CO.  The 
values  of  collisional  energy  exchange  coefficients  are  listed  in  Table  2. 

PHYSICS  OF  THF  PROGRAM  EXCHNG 


Aurora  and  Drizzle  Fluxes 

When  energetic  electrons  enter  into  the  atmosphere,  they  cause 
Ionization  and  initiate  many  excitation  proce  oes.  In  all  thes-  pro¬ 
cesses,  the  transfer  of  energy  occurs  from  ore  particle  to  another  and 
from  .ne  form  into  the  other.  The  excitation  in  the  vibration-rotation 
mode  thus  ceased  makes  important  contributions  to  the  infrared  emissions 
in  the  upper  atmosphere.  The  more  important  sources  of  these  energetic 
electrons  are  aurora  and  electron  drizzles. 

Three  "standards"  were  chosen  for  the  primary  electron  fluxes  and 
represent  IBC-I,  -II,  and  -III  auroras.  These  are  all  expressed  in  ex¬ 
ponential  form  as 

<KE)  =  0o-e"E/Eo  (5) 

where  E  =  the  energy  of  the  primary  electron, 

^o  *  the  l°w-energy  flux  limit  as  E  approaches  zero,  and 

Ec  =  t*1e  e-f' elding  energy. 

Table  3  lists  the  values  of  and  Eq  used  in  this  report  for  the  three 
IBC  auroras . 

Excepting  auroras,  there  is  an  additional  source  of  energetic  elec¬ 
trons  for  the  high-latitude  ionosphere,  which  we  term  as  an  electron 
drizzle.  The  satellite  measurement  made  by  Burch  [1968]  gave  evidence 
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TABLE  3 

THE  VALUES  OF  <p  AND  E  FOR  IBC-I 
o  o  9 


-II,  and  -III  AURORAS 


(electrons  cm-2  ster-1  sec-1  kev-1) 


Eq  (kev) 


I 

II 

III 


107 

108 

108 


3.5 


8.5 

12.5 


The  vaines  of  <f>?  and  E0  for  IBC-I  and  -II  were  given  by  Ulwiok  e '  al 

s^rr,  iz  SEn_ir  £■- ■*i  W  « 

Figure  19.  C  If  11 •  and  -m  auroras  are  plotted  in 


to  the  existence  of  an  electron  drizzle  and  indicated  that  its  energy 
spectrum  is  softer  than  that  of  most  auroras  and  ranged  approximately 
fi.om  50  ev  to  10  kev.  This  long  term  bombardment  of  energetic  electrons 
can  make  some  significant  contribution  to  the  ionization  and  excitation 
in  the  normal  high-latitude  ionosnhere  in  the  nighttime.  According  to 
BUr°h'8  measurement,  the  drizzle  flux  can  be  approximated  by 

<ME)  =  4.0  x  107  (E)-2  (electrons  cm-2  ster-1  sec-1  kev-1)  (6) 
with  E,  the  electron  energy,  ranging  from  0.05  kev  to  10  kev. 

Chemical  Model 

After  the  ionization  of  the  major  species  N2,  02  and  0  in  the 
ionosphere  by  the  primary  electron  flux,  the  ions  and  electrons  thus 
produced  will  immediately  undergo  a  series  of  secondary  reactions. 

The  primary  and  secondary  reactions  involved  in  variations  of  ion  and 
electron  concentrations  listed  below  are  the  important  ones  considered 
in  our  calculation. 

n2  +  (energetic  electron)  •*  n£  +  e 

°2  +  (energetic  electron)  -*■  of  +  e 

0  +  ''energetic  electron)  -*■  o+  +  e 
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0+  +  N2  ■+  N0+  +  N 

1.0  x  10"12  •  (T/300)-1 

0+  +  02  -*■  02  +  0 

2.0  x  10-11  •  (T/300)-*5 

0+  +  NO  -v  NO+  +  0 

2.4  x  10"11 

N2  +  0  -*■  0+  +  N2 
1.0  x  10“12 

N2  +  0  -*■  N0+  +  N(2D) 

1.4  x  10"10 

N2  +  02  -  02  +  N2 
5.0  x  10-11 

N2  +  NO  -*■  N0+  +  N2 
5.0  x  10-10 

02  +  N2  -*■  N0+  +  NO 
5.0  x  10“16 

0  2  +  NO  -►  N0+  +  02 
6.3  x  10-10 

02  +  N  -*■  N0+  +  0 

1.8  x  10“10 

e  +  02  ■►0  +  0 

2.1  x  10“7  •  (T/300)-1 

e  +  N2  >  N  +  N(2D) 

2.8  x  10-7  •  (T/300)"-33 


(7) 

\_Dunkin  et  al.  1968] 

(8) 

[ Dunkin  et  al.  1968] 

(9) 

[Fehsenfeld  et  al.  1965] 

(10) 

T Ferguson  1967] 

(ID 

[Fehsenfeld  et  al.  1970] 

(12) 

[Dunkin  et  al.  1968] 

(13) 

[Goldan  et  at.  1966] 

(14) 

[Ferguson  1967] 

(15) 

[Fehsenfeld  et  al.  1970] 

(16) 

[Goldan  et  al.  1966] 

(17) 

[Biondi  1969] 


(18) 


[Biondi  196 9] 
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e  +  NO  -v  N(2D)  +  0 

4.1  x  10“7  •  (T/300)-3 

N(2D)  +  02  >  NO  +  0 


6.0  x  10~^2 


N  +  NO  N2  +  0 
2.2  x  10_1 1 


[ Biondi  1969] 


[Lin  and  Kaufman  1970] 


[Phillips  and  Schiff  1962] 


This  reaction  set  does  not  include  the  excitation  processes  due  to 
secondary  electrons.  These  processes  are  taken  into  account  separately. 
The  only  neutral  reactions  included  in  this  model  are  N(2D)  +  02  and 
N  +  NO.  The  reaction  of  N(2D)  +  02  is  the  major  source  for  NO.  The 
production  of  N(2D)  comes  mainly  from  dissociative  recombination  reac¬ 
tions  of  N*  and  NO+.  According  to  Micheles  [1970],  the  production  of 
N(  S  :  D  :  2P)  from  the  reaction  of  Nf  +  e  is  in  the  ratio  40  =  40  =  10. 
The  quantum  yield  of  production  of  N(2D)  from  the  reaction  N0+  +  e  is 
assumed  to  be  unity  in  the  calculations  presented  here. 

Time-dependent  calculations  are  applied  to  this  chemical  model  to 
determine  the  concentrations  of  the  electrons  and  various  ions  with  the 
aurora  and  electron  drizzle  input.  Production  rates  of  N+  0+,  and  0+ 
due  to  prlmarv  electron  flux  of  aurora  and  drizzle  must  be  evaluated 
before  one  can  make  any  time-dependent  calculation  of  concentrations. 

Calculation  of  Production  Rates  of  nt.  0 t.  and  0+  due  to 
Primary  Electron  Flux  of  Aurora  and  Drizzle 


The  ionization  rate  coeffi  ient  due  to  primary  electron  flux  is  a 
function  of  altitude  because  energetic  electrons  with  different  ener¬ 
gies  have  different  limits  on  their  penetration  depth  into  the  atmo¬ 
sphere.  The  maximum  penetration  depth  for  electrons  with  energy  E  kev 
can  be  calculated  by  using  the  following  formula  [Gerard  1970] 


4.6  x  10~6  •  (E)!-65  (gm/cm2) 
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This  penetration  depth  is  given  as  the  mass  of  air  penetrated  through 
in  a  unit  area  of  column.  The  actual  altitude  corresponding  to  a  ver¬ 
tical  penetration  depth  in  the  atmosphere  can  be  found  in  a  functional 
relationship,  X  =  X(h),  where  h  is  the  altitude  and  X  is  the  penetra¬ 
tion  depth.  This  functional  relationship  is  tabulated  in  Table  4. 
Therefore,  the  minimum  energy,  E^,  needed  for  electrons  to  penetrate 
down  to  a  particular  altitude  can  be  found  by  using  the  above  formula 
and  Table  4.  The  linear  ionization  rate  at  an  altitude  h  due  to  pri¬ 
mary  energetic  electrons  with  energy  E  kev  is  given  by  Lazarev  [l?87]  as 


q(h’E)  =  3Ta f  d(x)p 

m 


(cm-3  sec-*  kev-*) 


(23) 


where  E(kev)  >  E  ,  (kev) 
—  min  ' 


and 


with 


<(>  =  the  electron  flux 
p  =  mass  density  of  air  at  altitude  h 

At  -  average  energy  dissipation  of  primary  electrons  for  a 

single  ionization  to  produce  an  electron-ion  pair,  34  ev 
needed  for  air 

=  the  solid  angle  extended  by  a  cone  with  a  half  vertex 
angle  of  80° 

d(x)  =  4.2X  exp  (-X2  -  x)  +  0.48  exp(-17 .4X* • 37) ,  called  the 
dimensionless  energy 

X  =  X(h)/X  (E) 
m 


The  total  production  rate  of  ions  is  obtained  by  the  integration 


Q(h)  =  /"  q(h,E)dE 


min 


(cm-3  sec-*) 


(24) 


The  individual  production  rates  for  nJ,  of,  and  0+  from  the  ionization 
of  N2.  02,  and  0  by  primary  electron  flux  are  given  as 


Qi(h)  =  oTTiT  +  /  n.  +  on  Q(h)  (cm"3  sec_1)  (25) 

«2  N2  02  02  00 

where  the  subscript  i  represents  the  ith  species  N2,  02,  or  0;  n's  are 
the  concentrations;  and  o's  are  the  ionization  cross  sections.  Values 
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TABLE  4 

THE  PENETRATION  DEPTH  AND  DENSITY  OF  THE  ATMOSPHERE  [ Rees  1963] 

h 

(km) 

X 

(gm  cm-2) 

P 

(gm  cm 

"3> 

60 

2.34  -  1 

3.04  - 

7 

C4 

1.37  -  1 

1.89  - 

7 

68 

7.76  -  2 

1.15  - 

7 

72 

4.20  -  2 

6.71  - 

8 

76 

2.17  -  2 

3.71  - 

8 

80 

1.08  -  2 

1.94  - 

8 

84 

5.16  -  3 

9.59  - 

9 

88 

2.46  -  3 

4.56  - 

9 

92 

1.19  -  3 

2.17  - 

9 

96 

5.88  -  4 

9.93  - 

10 

100 

3.10  -  4 

4.48  - 

ia 

108 

1.04  -  4 

1.32  - 

10 

114 

5.18  -  5 

5.70  - 

11 

120 

2.75  -  5 

2.61  - 

11 

126 

1.65  -  5 

1.30  - 

11 

132 

1.08-5 

7.05  - 

12 

140 

6.90  -  6 

3.47  - 

12 

150 

4.40  -  6 

1.70  - 

12 

162 

2.96  -  6 

9.00  - 

13 

192 

1.41  -  6 

2.94  - 

13 

230 

7.05  -  7 

1.16  - 

13 

252 

5.02  -  7 

7.60  - 

14 

276 

3.59  -  7 

5.00  - 

14 

300 

2.65  -  7 

3.42  - 

14 

Nott:  Entries  in  the  Table  N  ±  n  are  to  be  read  as  N  X  10-n 
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calculated  for  QLy  Qq^,  and  Q0  are  sent  into  the  program  SIXBIT  [ Adams 
and  Meg'1. 1,  1970].  In  this  program  the  ionization  reactions  due  to  the 
energetic  electrons  in  aurora  and  drizzle  are  treated  in  the  same  way 
as  other  chemical  and  photochemical  reactions  for  the  time-dependent 
calculation  of  concentrations  of  ions  and  electrons. 


Time-Dependent  Calculations  of  the  Electron  and  Ion  Concentrations 

The  generalized  reaction  kinetics  program  SIXBIT  was  originally 
developed  for  the  time-dependent  calculation  of  concentrations  of  re¬ 
actants  and  products  in  a  chemical  model.  This  program  was  modified 
to  calculate  temporal  variations  of  electron  and  ion  concentrations 
with  energetic  electron  input.  In  addition  to  the  subroutine  added 
for  the  calculation  of  production  rates  of  N2,  02,  and  0+  due  to  pri¬ 
mary  electron  flux,  new  subroutines  for  the  computation  of  secondary 
electron  flux,  the  estimation  of  the  electron  temperature,  vibrational 
excitation  and  de-excitation  rates  of  N2  and  02 ,  and  electronic  exci¬ 
tation  rate  of  O^D)  were  also  added.  The  modified  program  is  called 
EXCHNG . 

All  the  following  calculations  are  made  for  altitude  t  nges  from 
60  km  to  150  km  in  the  nighttime.  Several  steps  are  taken  in  the  cal¬ 
culations.  The  electron  and  ion  concentrations  for  quiescent  condi¬ 
tions  in  the  nighttime  were  taken  from  Jones  and  Rees  [1972]  as  the 
initial  profiles  and  the  calculation  were  made  with  constant  drizzle 
input.  The  time-dependent  calculation  was  made  to  continue  until  quasi 
steady  state  for  the  concentrations  was  reached.  In  the  next  step,  the 
quasi  steady-state  profiles  with  drizzle  input  were  taken  as  the  initial 
profiles  in  time-dependent  calculations  with  constant  aurora  input.  The 
quasi  steady  state  is  also  reached  for  this  calculation.  Following  this 
procedure  the  results  obtained  for  IBC-II  and  -III  auroras  are  plotted 
in  Figures  20  and  21.  A  separate  calculation  was  made  with  constant 
IBC-II  aurora  input  starting  with  the  Initial  quiescent  profiles  from 
Jones  and  Rees.  The  initial  profiles  and  the  steady-state  profiles  with 
constant  aurora  input  are  plotted  in  Figure  22. 
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Calculation  of  Secondary  Electron  Flux 

This  calculation  follows  the  method  suggested  by  Rees  et  at.  [1969]. 
The  production  rate  of  secondary  electrons  is  evaluated  by  the  following 
expression 


n<Es) 


G(V 

2.53 


QN2  + 


C(Wn  ) 
2.77 


C(w  ) 

qo2  +  -E7T  (26) 

(electrons  cm  3  sec-1  ev_1) 


where  Q,  is  the  rate  of  ionization  due  to  primary  electrons,  and  E_  is 
1  s 

the  energy  of  the  secondary  electrons .  G(W)  is  a  function  of  W  given 
as 


G(W)  = 


ex 


p[- 


W 


31.5 


exp(- 


W 


with  W  =  I  +  E 


,  ,31. 6  +  A000  -  WN 

x  ln  ( -  - )  (27) 

31.6  -  A000  -  W 
(ev*1) 


and  I  =  the  ionization  potential,  15.6  ev  for  N2 ,  12.1  ev  for  02 ,  and 
13.6  ev  for  0. 


With  the  assumption  that  energy  loss  is  a  continuous  function  of  energy, 
the  differential  flux  of  secondary  electrons  is 


„CEs)  dEt 


*<ES>  -rm 


,  .  +  L(E) 

elect.  neut. 


(electrons  cm-2  sec"1  ev  1 )  (28) 


where  L(E)e^ect  and  L(E)  are  the  terms  of  energy  loss  to  ambient 

electrons  and  to  neutrals,  respectively,  per  unit  distance  traveled  by 
the  secondary  electrons.  The  function  L(E)  .  is  given  in  Rees' 

616C t • 

paper  as 


L(E) 


elect . 


1.95  x  IQ"'2  n(e) 


(ev  cm-1) 


(29) 


where  n(e)  is  the  ambient  electron  density  in  units  of  electrons/cm3. 

For  the  function  L(E)  ,  the  approximate  analytical  expression  for 

neut  • 

E  down  to  25  ev  is  given  by  Green  and  Peterson  [1968]  as 


L(E) 


neut. 


2T  "  t(f>  +  <f>  +  <!>  T  ■"! 


(ev  cm-1) 


(30) 
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where  qQ  =  6.514  x  10“11  ev2  cm2 

R  =  13.6  ev 
e 

and  0,  M,  A,  A,  B,  C,  and  g  are  the  parameters  which  are  listed  in 
Table  5  for  N2 ,  02,  and  0,  and  nA  is  the  concentration  of  ith  species, 
N2,  02,  or  0.  For  E  in  the  range  from  25  ev  to  1  ev,  Green's  approxi¬ 
mation  [1972]  is  used  in  the  calculation  of  L(E)  .  The  energy  of 

the  secondary  electrons  in  this  range  is  mainly  lost  to  electronic  and 
vibrational  excitation,  the  Rydberg  series,  and  the  ionization  continuum 
of  N2,  02,  and  0  in  the  upper  atmosphere.  Therefore,  the  energy  loss 
functions  will  be  discussed  individually  in  these  three  categories. 

TABLE  5 

PARAMETERS  USED  IN  GREEN  AND  PETERSON'S  APPROXIMATE  ANALYTIC 
EXPRESSION  OF  THE  ELECTRON  ENERGY  LOSS  TO  N2 ,  02  ,  AND  0 
[ Green  and  Feierson  1968] 


Parameter 

N2 

°2 

0 

ft 

0.783 

0.771 

0.758 

M 

-0.328 

-0.518 

-0.500 

A 

-3.14 

-3.59 

-3.5 

A  (ev) 

100 

100 

100 

B  (ev) 

217 

123 

765 

C  (ev) 

54.6 

54.1 

50 

_+ 

Z 

14.1 

15.2 

7.15 

1.  Loss  due 

to  electronic  and  vibrational 

excitation 

L(E)  = 

e 

l  .o.,(E)  for  E  >  W , 

ij  1  —  ij 

(ev  cm-1) 

(31) 

where  E  =  the  energy  variable,  which  is  the  energy  of  the  secondary 
electrons  in  the  discussion  considered  here 

subscript  i  =  the  ith  species,  N2,  02,  or  0 

subscript  j  =  the  jth  vibrational  level  or  electronic  state 

=  the  concentration  of  the  ith  species 


A3 


W  =  the  excitation  energy  of  the  jth  vibrational  level  or 
electronic  state  of  the  ith  species  relative  to  the 
lowest  vibrational  level  in  the  ground  electronic  state 

o..(E)  =  the  excitation  cross  section  of  the  jth  vibrational 
level  or  electronic  state  of  the  ith  species 

The  excitation  cross  section  has  an  expression 


q  (f  C  )  W, .  .  W. .  0 

°ij(E)  -  t-?11  [1  - 


(cm2) 


ij 


(32) 


where  q  =  4ira  2Ry2  =  6.51A  x  10” * 4  cm2  ev2 
o  o 

(aQ  is  the  Bohr  radius  and  Ry  is  the  Rydberg  energy,  13.6  ev)  , 
and  fQCo,  6,  and  v  are  the  input  parameters  which  have  differ¬ 

ent  values  for  the  various  vibrational  levels  or  electronic  states 
for  different  species.  All  parameters  for  the  different  vibrational 
levels  or  electronic  states  of  N2  >  O2 ,  and  0  are  listed  in  Table  6. 


2.  Loss  due  to  Rydberg  series 

L(E)R‘EljniWViJ(E> 


(ev  cm'  1 ) 


(33) 


where  E,  the  subscripts  1  and  j,  n^  and  (E)  have  the  same  meaning 
as  those  in  category  (1)  but  WI^  is  the  ionization  potential  for  the 
jth  Rydberg  series  of  the  ith  species. 

The  expression  for  the  cross  section  has  the  form 


^nf-1)  (f  C  ) 


W, 


W. 


°ij  (E) 


£ 

n=n , 


W 


ijn 


o_e  [j  .  (-iis,®]'1 


ijR 


(cm2) 


(34) 


where  q  =  6.514  x  10” 14  cm2  ev2 
o 


^o^o^n  (n  -  6)  3 


(ni  £  n  <  nf) 


(foCo*R  =  2(nf  -6  -  h)z 
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tA5l£  6 

THE  PARAMETERS  USED  IN  CALCULATION  OF  THE  ENERl  Y  LOSS  OF 
SECONDARY  ELECTRONS  TO  NEUTRALS  FUF  TO  VIBRATIONAL  AND  ELECTRONIC 

EXCITATION  [Green  1972] 


Excited  State 

Wjj (ev) 

6 

V 

W.  j  (e  v) 

0 

f  C 
o  o 

n2 

N2(v=1-3) 

1.850 

1.000 

1.000 

0.570 

7.000 

0.273200 

N2 (v=4-8) 

2.150 

1.000 

1.000 

1  .680 

9.000 

0.241300 

Ah+ 

u 

6.168 

1.000 

1.000 

6.168 

3.000 

0.226000 

B  3ir 

g 

7.353 

3.000 

1.000 

7.353 

3.000 

0.177500 

C3Tru 

11.031 

3.000 

1.000 

11.031 

3.000 

0.281200 

E’E8 

11.900 

3.000 

1  .000 

11.900 

3.000 

0.047800 

8.548 

1.000 

1.000 

8.548 

1.000 

0.136200 

a‘Eg 

12.250 

2.300 

1 .000 

12.250 

1.000 

0.027200 

12.500 

1.000 

3.000 

12.800 

0.750 

0.672000 

b1^ 

13.300 

1.000 

3.000 

14.000 

0.750 

0.334000 

2p2(3P)3s  4P 

21.000 

1.000 

2.000 

24.000 

0.750 

0.200000 

2p2(]D)3s  3P 

21.000 

1 .000 

2.000 

24.000 

0.750 

0.047000 

2p2 ( 3P) 3s  2P 

21.000 

1.000 

2.000 

23.000 

0.800 

0.085000 

2p2(1D)3s  2D 

22.000 

1.000 

2.000 

25.000 

0.850 

0.064000 

02 

9.900 

1.000 

3.000 

9.900 

0.850 

0.080000 

alAg 

0.980 

1.000 

3.000 

0.980 

0.750 

0.000990 

b'Eg 

1.640 

3.000 

1.000 

1.640 

3.000 

0.000500 

4.500 

1.000 

1.000 

4.500 

0.900 

0.021000 

8.400 

1.000 

2.000 

8.400 

0.750 

0.230000 

16.100 

1.000 

3.000 

19.000 

0.900 

0.067000 

19.200 

1.000 

3.000 

22.000 

0.750 

0.016000 

17.600 

1.000 

3.000 

21.000 

1.000 

0.032000 

14.600 

1.000 

3.000 

18.000 

0.750 

0.085000 
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TABLE  6  (cont.) 


Excited  State 

Wy  (ev) 

8 

V 

Wij (ev) 

n 

f  c 

o  o 

0 

0(!d) 

0(JS) 

1.850 

4.180 

1.000 

0.500 

2.000 

1.000 

1.960 

4.180 

1.000 

1.000 

0.010000 

0.004200 

Wijn  "  WIij  '  (n  -y«5)^(ev) 

(Ry  =  Rydberg  energv) 

Wij»"  <UlJnf  +WI1J>/2 

* 

and  fQ,  <5,  ni,  nf,  fi,  8,  and  v  are  the  input  parameters.  In  the  above 
expression  of  (E) ,  the  last  term  represents  an  effective  state  which 

is  the  result  of  the  summation  over  all  the  terms  with  n  >_  nf .  The 
parameters  needed  in  the  calculation  of  this  loss  function  are  given  in 
Table  7.  nj_  and  nf  have  the  values  3  and  5  respectively  for  the  above. 

TABLE  7 

THE  PARAMETERS  USED  IN  CALCULATION  OF  ENERGY  LOSS  OF  SECONDARY 
ELECTRONS  TO  NEUTRALS  DUE  TO  THE  EXCITATION  OF  RYDBERG  SERIES 

[ Green  1972] 


Rydberg  Series 

WIjj  (ev) 

6 

V 

6 

0 

f* 

0 

n2 

x24 

15.580 

1.000 

3.000 

0.700 

0.750 

14.349999 

A27Tu 

16.730 

1.000 

3.000 

1.040 

0.750 

1.646000 

B2^ 

18.750 

1.000 

3.000 

0.870 

0.750 

1.000000 

D2TTg 

22.000 

1.000 

3.000 

1.530 

0.750 

0 . 700000 

C2I+ 

23.600 

1.000 

3.000 

0.800 

0.750 

0.700000 
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TABLE  7  (cont.) 


Rydberg  Series 

Wly  (ev) 

S 

V 

6 

ft 

f* 

0 

°2 

x2’8 

12.100 

1.000 

3.000 

1.080 

0.750 

1 . 200000 

a4"u 

16.100 

1.000 

3.000 

1.000 

0.750 

2.500000 

A2”  u 

16.900 

1.000 

3.000 

1.070 

0.750 

2.099999 

b"E- 

18 . 200 

1.000 

3.000 

0.690 

0.750 

1.799999 

n 

23.000 

1.000 

3.000 

1.000 

0.750 

0 . 700000 

u 

13.600 

1.000 

3.000 

1.160 

0.750 

0.289500 

13.600 

1.000 

2.000 

1.240 

3.000 

0.270200 

13.600 

1.000 

1.000 

0.810 

3.000 

0.077200 

13.600 

2.000 

1.000 

0.690 

0.750 

0.011600 

13.600 

1.000 

3.000 

0.010 

0.750 

0.400500 

13.600 

1.000 

2.000 

0.010 

3.000 

0.373800 

16.900 

1.000 

3.000 

1.210 

0.750 

0.481500 

16.900 

1.000 

2.000 

1.180 

3.000 

0.449400 

16.900 

2.000 

1.000 

0.920 

0.750 

0.019300 

16.900 

1.000 

1.000 

0.920 

3.000 

0.128400 

16.900 

2.000 

1.000 

0.810 

0.750 

0.019300 

16.900 

]  .000 

1.000 

0.780 

3.000 

0.128400 

16.900 

2.000 

1.000 

0.780 

0.750 

0.019300 

16.900 

1.000 

1.000 

0.780 

3.000 

0.128400 

16.900 

1.000 

3  000 

0.040 

0.750 

0.298500 

16.900 

1.000 

2.000 

0.040 

3.000 

0.278600 

16.900 

1.000 

3.000 

0.040 

0.750 

0.459000 

16.900 

1.000 

2.000 

0.040 

3.000 

0.428400 

16.900 

1.000 

3.000 

0.040 

0.750 

0.237000 

16.900 

1.000 

2.000 

0.040 

3.000 

0.221200 

18.500 

1.000 

3.000 

1.250 

0.750 

0.297000 

18.500 

1.000 

2.000 

1.190 

3.000 

0.277200 

18.500 

2.000 

1.000 

1.040 

0.750 

0.011900 

18.500 

1.000 

1.000 

1.040 

3.000 

0.079200 
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TABLE  7  (cont.) 


Rydberg  Series  Wl^ev) 

8 

V 

<S 

0 

f* 

o 

18.500 

2.000 

1.000 

0.770 

0.750 

0.011900 

18.500 

1.000 

1.000 

0.770 

3.000 

0.079200 

18.500 

2.000 

1.000 

0.770 

0.750 

0  011900 

18.500 

1.000 

1.000 

0.730 

3.000 

0.079200 

18.500 

1.000 

3.000 

0.050 

0.750 

0.235500 

18.500 

1.000 

2.000 

0.050 

3.000 

0.219800 

18.500 

1.000 

3.000 

0.050 

0.750 

0.678000 

18.500 

1.000 

2.000 

0.050 

3.000 

0.632800 

3.  Loss  due  to  ionization  cont i nuum 


HE). 


In J 
ij  10 


M 


W^S^CE.^dT 


(ev  cm-1) 


(35) 


where  E,  the  subscripts  i  and  j,  n^  have  the  same  meaning  as  those  in 

categories  (1)  and  (2),  but  «  1^  +  T  with  1^  representing  the 

ionization  energy  for  the  jth  continuum  state  of  the  ith  species. 

T  =  kinetic  energy  of  the  electron  released  from  the  ionized  mole¬ 
cule 

(E  -  the  possible  maximum  value  of  T 

Sij(E,T)  =  the  differential  ionization  cross  section  for  a  continuum 
state  j  of  the  ith  species 


The  integration  in  the  above  expression  for  L(E)^  yields 


M  I  a  E  +  I 

S^CE.TMT  =  qoAo(-|i)P  if  [(— - 

?  r 


W 


ij  ij 


2E 


)  *r  -  (-11)  V] 


(ev  cm2) 

(36) 


where  q 

o 


a 


r 


6.514  x  1CT14  cm2  ev2 
v  r 

cr (~1)  with  Cr  as  the  number  of  combinations  of  v  distinct 
things  taken  r  at  a  time 


♦  *  fl  -  P  +  T 
r  r 

Tr  Br 


48 


an<^  an<^  v  are  the  input  parameters. 

Table  8  presents  the  values  of  various  parameters  needed  in  the  calcu¬ 
lation  of  L(E)i. 

The  total  energy  loss  of  the  secondary  electrons  to  neutrals  within 

the  energy  range  from  25  ev  to  1  ev  is  the  summation  of  the  above  three 
losses . 

L^neut.  =  L^e  +  L^R  +  L(E>i  (ev  cm-1)  (37) 


TABLE  8 

PARAMETERS  USED  IN  THE  CALCULATION  OF  ENERGY  LOSS  OF  SECONDARY 
ELECTRONS  TO  NEUTRALS  DUE  TO  THE  IONIZATION  CONTINUA  [Green  1972] 


Ion  State 

rij 

8 

V 

P 

0 

A  ( ev  1 ) 
o  ' 

n2 

*3 

15.58 

1.00 

2.0 

1.20 

0.80 

0.37000 

A2iru 

16.73 

1.00 

1.0 

1.20 

0.83 

0.16000 

18.75 

0.15 

1.0 

1.20 

0.94 

0.36800 

d2’8 

22.00 

1.00 

2.0 

1.20 

0.83 

0.05600 

c2ij 

23.60 

1.00 

2.0 

1.20 

0.83 

0.06000 

o2 

x\ 

12.10 

1.00 

2.0 

1.10 

0.80 

0.05800 

a4TTu 

16.10 

1.00 

2.0 

1.10 

0.80 

0.15000 

A2ttu 

16.90 

1.00 

2.0 

1.10 

0.80 

0.15000 

b“i- 

18.20 

1.00 

2.0 

1.10 

0.80 

0.13000 

23.00 

1.00 

2.0 

1.10 

0.80 

0.06400 

18.00 

1.00 

3.0 

1.10 

0.93 

0.40000 

22.00 

1.00 

3.0 

1.10 

0.93 

0.25000 

0 

13.60 

1.00 

2.0 

1.20 

0.85 

0.26500 

16.90 

1.00 

2.0 

1.20 

0.85 

0.33000 

18.50 

1.00 

2.0 

1.20 

0.85 

0.18000 
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gj  the  Electron  Temperature,  Vibrational  Excitation,  and 
De-excitation  Rates  of  N?  and  0?  and  the  Electronic  Excitation  Rate  of 

0(1D) 

After  the  secondary  electron  flux  is  determined,  the  electron 
temperature  can  be  calculated  by  making  the  energy  input  to  the  am¬ 
bient  electrons  by  the  secondary  electrons 

00 

Ein  =  'rc)L^Es^elec.  dEs  (ev  cm~3  sec_1)  (38) 

equal  to  the  energy  loss  from  the  heated  ambient  electrons  to  the  neutral 
molecules.  This  energy  loss  is  calculated  by  following  the  method 
suggested  by  Megill  [1965]. 

Eloss  =  v(u)  8(u) (u  -  uo)  [M][E]  (ev  cm"  3  sec"1)  (39) 

where  u  =  average  energy  of  the  ambient  electrons,  3/2kT  ,  with  Ie  as 
the  electron  temperature  and  k  as  the  Boltzmaneconstant 

uo  =  ayerage  thermal  energy  of  the  neutrals,  3/2kT„  with  Tn  as 
the  temperature  of  neutral  particles 

g(u)  =  average  fractional  energy  loss  per  collision 

v(u)  -  average  collision  frequency  between  ambient  electrons  and 
neutral  particles  normalized  to  a  neutral  particle  density 
of  one  per  cmJ 

[m]  -  total  number  density  of  neutral  particles 
[E]  =  number  density  of  the  ambient  electrons 
Both  g(u)  and  v(u)  are  functions  of  the  average  energy  of  ambient  elec¬ 
trons  and  are  plotted  in  Figure  23.  The  condition  of  steady  state  is 
assumed  for  the  average  energy  of  ambient  electrons,  that  is 

Ein  =  Eloss 

The  electron  temperature,  Te,  which  is  included  in  the  above  equation 
through  u,  can  be  solved  numerically  from  this  equation. 

Calculation  of  the  excitation  rates  of  N2  and  02  into  the  first 
vibrational  level  in  their  ground  state  by  collision  with  secondary 
electrons  is  simply  made  by  the  integrations 

Vex(N2)  =  ^WV  dEs  (sec-1)  (40) 


lal  energy  loss  g  plotted 
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Vex^°2^  =  *°02 (Es>  dEs  (sec-1)  (41) 

where  o^^)  is  the  excitation  cross  section  for  N2  given  by  Chen 
[1964]  end  °q2 (Eg)  is  the  one  for  02  given  by  Linder  and  Sahmidt  [l97l]. 
A  similar  integration  yields  the  excitation  rate  of  0  into  its 
state. 

Veji<0(>D))  -  /*(Es).„0(1[))(Es)  dEs  (sec-1)  (42) 


where  a0(lD)  is  the  excitation  cross  section  for  0(JD). 

An  expression  for  this  cross  section  is  given  by  Green  [1972]  as 


C70(1D)(E')  = 


q  (f  C  ) 

— w-2-°  (V  [i.o  -  (yr 


(cm2) 


(43) 


where  qQ  =  6.514  x  10~14  cm2  ev2 
f  C  =  0.01 

o  o 

w  =  1.85  ev 
0-1.00 
B  =  1.00 
v  -  2.00 

The  deactivation  rate  for  the  vibrational  level  (v  =  1)  of  N2  by 
collisions  with  slow  ambient  electrons,  Kd>  is  evaluated  by  using  the 
analytic  approximation  to  the  cross  section  given  by  Chen  [1964],  as 
shown  below  for  different  energy  intervals 


O]  5  x  10  13  (1  -  ^1)  exp  (-  ~)  (cm2)  (1.5  ev  <_  E) 


(44) 


o2  -  0.3  x  10  16  exp  (-  ~)  (cm2)  (0  _<  E  0.2  ev) 


(45) 


o3  -  0.05  x  10-16  (cm2)  (0.2  ev  <  E  <  1.5  ev)  (46 

where  o’s  are  given  in  the  units  of  cm2  and  E  is  the  kinetic  energy  of 
ambient  electrons. 

The  sum  over  the  three  energy  intervals  gives 


o  =  O]  +  a2  +  a3  (cm2) 
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The  deacti nation  rate  is  given  by 

=  /va(E)  f(E)  dE  (cm3  sec-3)  (47) 

where  f(E)  =  (^)**  (kT,)"  3/* .v.e“E/kTe  (ev-3) 

(the  Maxwell  distribution  of  kinetic  energy) 
v  =  /2E/rae  (the  velocity  of  ambient  electrons) 

me  =  electron  mass 

Tg  =  electron  temperature 
k  =  Boltman  constant 

Finally,  the  above  integration  gives  rise  to  the  following  expression 
for  Kd 

(T  )** 

Kd  =  3,34  X  10  5  711 59  +  T  )2  +  3,11  x  10“12(T  )**  (cm3  sec*3)  (48) 

e  e 

ELECTRIC  FIELD  EXCITATION  OF  INFRARED  EMISSIONS 

Recent  measurements  [Westcott  et  al.  1969,  Hepner  et  al.  1971]  indi¬ 
cate  that  the  auroral  regions  are  quite  often  subjected  to  the  presence 
of  relatively  large  electric  fields.  These  fields  are  generally  meridi¬ 
onal  in  nature  and  vary  considerably  with  time.  Fields  as  high  as  100  mv/ 
meter  are  fairly  often  observed.  A  field  as  large  as  200  mv/meter  would 
be  very  large  when  a  field  of  10  mv/meter  is  near  what  one  would  expect 
for  "quiet"  times. 

The  geographic  extent  of  these  fields  is  at  present  uncertain,  ex¬ 
cept  that  at  times  fields  have  been  observed  at  widely  separate  points; 
e.g..  Point  Barrow,  Alaska,  and  Fairbanks,  Alaska.  One  could,  therefore, 
expect  that  these  field  structures  are  of  great  extent  compared  to  auro¬ 
ral  features  —  hundreds  of  kilometers  as  compared  to  several  kilometers 
for  auroral  forms.  It  is  not  clear  at  present  whether  or  not  these 
fields  exhibit  the  characteristics  of  current  sources  or  voltage  sources. 
In  other  words,  does  the  electric  field  remain  constant  in  an  auroral 
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arc  when  the  ionization  and  consequently  the  conductivity  is  increased 
greatly?  There  are  no  definitive  answers  to  this  question  at  present, 
although  relatively  simple  arguments  concerning  the  total  amount  of 
charge  available  in  the  magnetosphere  would  argue  for  a  current  source. 

It  is  of  interest  to  consider  the  energy  inputs  of  a  system  such 
as  this.  Heaps  [1972]  and  Harris  [1972]  have  done  this  by  studying 
the  ohmic  heat  input  to  the  atmosphere  as  well  the  ion  drag  and 
some  of  the  motions  consequent  to  these  energy  and  momentum  inputs. 

In  this  section  we  explore  the  way  in  which  energy  is  placed  into 
vibrationally  excited  states  of  the  02  and  N2  molecules  by  collisions 
with  ions  driven  by  the  electric  fields.  We  have  assumed  that  the 
process 


A+  + 


BC 


A+  + 


BC 


(49) 


is  the  only  type  of  process  that  proceeds.  We  know  that  the  process 
is  much  more  complicated  than  this,  because  in  drifting  experiments  in 
which  ions  drift  in  their  parent  gases,  much  more  complicated  processes 
occur.  For  example,  N2  drifting  in  N2  results  in  the  ions  N+,  N2,  N+, 
and  N4  [ Bloomfield  and  Hasted  1966].  One  item  that  is  quite  clear  is 
that  processes  like  this  are  continually  breaking  and  making  bonds. 
These  are  processes  which  are  likely  to  create  significant  amounts  of 
vibrational  energy.  To  the  knowledge  of  these  authors,  no  quantita¬ 
tive  information  exists  as  to  cross  sections  for  processes  like  equa¬ 
tion  (49).  We  have,  therefore,  adopted  a  set  of  cross  sections  which 
we  think  are  possible.  Indeed,  we  feel  that  these  are  on  the  conserva¬ 
tive  side  in  the  sense  that  our  estimates  are  probably  smaller  than 
the  cross  sections  themselves .  We  have  considered  excitation  of 
02(v=l),  N2 (v=1) ,  and  02(1A  )  by  0+,  nJ,  N0+,  and  02  ions.  The  cross 

o 

sections  assumed  for  each  of  the  processes  is  given  in  Table  9. 

The  calculation  that  we  have  made  is  considered  to  be  a  very  con¬ 
servative  one.  Conservative  again  in  the  sense  that  we  expect  that  we 
calculate  excitation  rates  smaller  than  those  which  may  actually  occur. 
We  assume  that  the  ion  is  accelerated  by  the  electric  field  until  it 
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TABLE  9 

CROSS  SECTIONS  IN  cm2 


STATE  EXCITED 


collides  with  another  molecule  or  until  it  is  forced  into  cyclical 
motion  by  the  magnetic  field.  The  calculations  assume  a  mean  free 
path  determined  by 

L  ‘  [5F  (50) 

where  [M]  Is  the  sum  of  the  0,  02 ,  and  »2  densities.  The  cross  section 
is  assumed  to  be  10-14  for  all  species  and  ions  for  purposes  of  mean 
free  path  calculations.  This  cross  section  is  analogous  to  a  diffusion 
cross  section  for  electrons.  If  the  mean  free  path  is  longer  than  a 
gyro-radius  of  the  ion,  then  the  maximum  energy  which  the  ion  can  at¬ 
tain  is  determined  by  the  magnetic  Held.  At  each  altitude,  both  the 
mean  free  path  and  the  gyro-radius  are  calculated  for  each  ion.  The 
energy  as  a  function  of  time  is  then  calculated  using  the  appropriate 
equations.  In  all  cases  the  ions  are  considered  to  start  from  a  zero 
velocity.  The  fraction  of  cheir  energy  that  is  spent  above  the  thresh¬ 
old  energy  of  the  state  berlng  considered  is  then  calculated.  The  ap¬ 
propriate  equations  are  given  in  Table  10.  For  example,  in  considering 
the  excitation  of  N2(v=l),  we  use  a  threshold  of  .3  ev.  The  average 

velocity  is  also  calculated.  V7e  then  form  the  rate  of  excitation  per 
molecul®  as 

R(j)  = 


v4[n  ]  T  (i,j) 
_i  -if 


(51) 


TABLE  10 

ION  ENERGY  CAI CULATIONS 
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Physical  Paramete-s 


Equations 


Gyro-radius  of  ion  r  =  eE/Mw2 

Mean  free  path  of  ion  d  =  l/(o*[lon]) 


Average  velocity  of  ions 

d  >  R 

d  <  R 

AVEL  =  1.414E/B 

AVEL 

-  */uED/2M 

Average  energy  of  ion 

AENR  =  ME2/B2 

AENR 

=  eEd/4 

Maximum  energy  of  ion 

MENR  =  2ME2 /B2 

MENR 

=  eEd 

Fractional  time  ion 
have  energy  >  EXR 

TA  =  (l/7r)sin-1 
(K.E./ME2/B2)-1) 

TA 

=  1  -  ED /MENR 

Excitation  rate  per 
molecule 

EXR  (0  •  [  Ion]  TA 

AVEL) 

for  each  ion 

where  w  =  eB/M 

e  =  electronic  charge 
E  =  electric  field  strength 
B  =  magnetic  field  strength 
[ion]  =  concentration  of  ions 
a  cross  section 
M  £  mass  of  ion 


where  j  represents  the  state  being  considered;  a(i,j)  is  obtained  from 
Table  9;  vi  and  [n^  are  the  average  velocity  and  number  density  of  the 
ith  ion;  Tf(i,j)  is  the  fraction  of  its  time  which  the  ith  ion  spends 
with  enough  energy  to  excite  the  jth  state. 

The  results  of  these  calculations  are  then  entered  as  excitation 
ra.es  of  the  N2  and  O2  vib  ational  states  in  the  program  previously 
described,  and  the  resulting  infrared  emissions  are  calculated.  The 
O2 ( excitation  rate  was  also  calculated  using  this  technique,  but 
the  results  will  not  be  presented  in  this  report.  Figures  24,  25,  and 
26  show  the  results  of  these  calculations  for  electric  fields  of  0.05, 


intensity  of  0.1  volts .meter 


59 


0.1,  and  0.2  volts/meter.  By  referring  to  Figure  77,  one  can  see  the 
marked  affect  of  the  electric  field  on  N2  vibrational  temperature. 

This  is  a  plot  of  the  steady  state  ratio  of  N2(v=l)  as  opposed  to  that 

which  would  be  calculated  from  the  assumption  of  a  Boltzmann  distribu¬ 
tion. 

As  was  stated  earlier,  a  number  of  conservative  assumptions  have 
been  made  in  these  calculations.  We  have  ignored  the  initial  ion  ener¬ 
gies  which,  for  auroral  altitudes,  may  represent  a  maxwellian  distribu¬ 
tion  function  with  average  energies  of  several  hundredths  of  an  elec¬ 
tron  volt.  When  one  is  considering  the  excitation  of  states  as  low  as 
0.2  ev  (for  excitation  of  02(v=l)),  inclusion  of  this  energy  would  in¬ 
crease  the  calculated  excitation  rate  substantially.  A  larger  term 
that  has  been  ignored  is  the  substantial  increase  in  random  kinetic 
energy  as  measured  in  the  drift  frame  of  reference.  This  energy  can, 
depending  on  the  relative  mass  of  the  ions,  become  substantial  if  only 
elastic  collisions  occur  {Megill  and  Carle  ten  1964].  Inelastic  pro¬ 
cesses  may  well  dominate  in  the  auroral  situation  so  that  accurate 
evaluation  of  this  term  is  difficult.  The  effect  of  neglecting  this 
term  is,  again,  to  underestimate  the  intensities. 

IR  BAND  MODEL  PROGRAM  (PROGRAM  BCKGND)  PHYSICS 

A  very  general  description  of  the  methods  used  in  computing  the 
band  radiance  for  each  of  tne  infrared  bands  under  consideration  was 
given  in  a  previous  section  of  this  report.  It  is  the  purpose  of  this 
section  to  consider  the  physics  of  this  program  in  more  depth. 

-ince  this  research  used  the  IR  band  model  program  of  Corbin  et  al. 
[1970]  as  a  starting  point,  much  of  the  physics  described  in  their  re¬ 
port  will  also  be  applicable  here.  On  the  other  hand,  many  new  ideas 
and  the  associated  physical  processes  have  been  incorporated  in  the 
original  work  with  the  consequent  change  in  overall  philosophy.  There¬ 
fore,  a  complete  physical  description  will  be  made  even  though  it  may 
mean  duplicating  some  of  the  physical  description  found  in  the  report 
by  Corbin  et  al.  [1970]. 
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General  Concepts 

Ultimately,  it  is  the  population  of  the  various  vibrational  states 
of  a  molecule  which  determines  the  infrared  radiation  from  that  mole¬ 
cule.  Molecules  may  be  vibrationally  excited  by  one  of  several  means; 
absorption  of  radiation,  fluorescent  excitation,  chemiluminescent  pro¬ 
cesses,  or  by  one  of  several  collisional  excitation  processes.  These 
processes,  in  any  kind  of  steady  state,  are  balanced  by  radiative  or 
collisional  de-excitation. 

In  the  lower  regions  of  the  atmosphere  where  collisional  frequen¬ 
cies  are  high,  collisional  excitation  and  de-excitation  processes  are 
the  more  important.  As  one  proceeds  higher  in  the  atmosphere,  colli¬ 
sional  frequencies  decrease  and  radiative  processes  become  more  and 
more  important.  In  order  for  a  model  to  predict  IR  band  radiance  over 
a  wide  altitude  region,  it  is  necessary  to  include  all  processes  over 
all  altitudes.  This  is  especially  true  when  one  considers  the  varia¬ 
tion  of  excitation  and  de-excitation  cross-sections  for  the  same  pro¬ 
cesses  from  species  to  species. 

In  the  proceeding  sections  it  was  shown  that  electrons  and/or 
electric  fields  tend  to  put  a  substantial  amount  of  their  energy  into 
the  vibrational  energy  of  N2  and  02  as  well  as  into  0(*D).  Physically 
this  energy  input  has  bean  computed  in  the  form  of  production  rates. 
Additionally,  the  reverse  processes  also  take  place  for  N2  and  02 . 

That  is,  a  certain  amount  of  vibrational  energy  is  lost  to  the  ambi¬ 
ent  electrons.  For  these  processes,  the  computational  input  is  in  the 
form  of  a  loss  rate.  In  the  case  of  O^D),  much  of  this  energy  is 
transferred  to  N2  and  some,  of  course,  to  02. 

These  production  and  loss  rates  for  N2(v=l),  02(v=l),  and  O^D) 
are  obtained  from  the  computer  program  EXCHNG  d  scribed  earlier. 

The  values  obtained  depend  upon  the  auroral  classification  used.  A 
primary  electron  flux  was  assumed  as  representative  of  an  IBC-I,  -II, 
and  -III  aurora  and  from  this  flux  the  production  and  loss  rates  were 
determined.  The  auroral  classifications  and  representative  primary 
electron  flux  are  listed  in  a  preceding  section.  In  addition  to  these 
aurorally  disturbed  conditions,  three  categories  of  "quiescent”  condi- 
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tions  were  examined.  These  categories  included  a  totally  quiet  condi¬ 
tion,  a  condition  whereby  an  electron  drizzle  of  flux  4  x  107  x  E-2 
cm  2  sec-1  ster-1  kev-1  exists,  and  a  condition  which  is  totally  quiet 
but  in  which  an  electric  field  exists  in  the  ionosphere.  Field  values 
of  0.05,  0.1,  and  0.2  volts/meter  were  used  as  representative.  The 
electric  fields,  as  in  the  case  of  electron  precipitation,  tend  to  put 
energy  into  the  vibrational  states  of  N2  and  02  since  these  are  the 
major  constituents.  Production  rates  for  N2(v=l)  and  02(v=l)  resulting 

from  the  electric  fields  were  computed  as  described  in  a  previous  sec¬ 
tion. 

Atmospheric  disturbances  represent  additional  energy  inputs  to 
the  atmospheric  system.  It  is  the  fact  that  a  considerable  amount  of 
this  energy  ends  up  as  vibrational  energy  in  N2  and  02  that  makes  a 
close  examination  of  the  energy  exchange  processes  of  these  vibration- 
ally  excited  molecules  important.  Further,  it  makes  the  near  resonant 
V-V  exchange  process  probably  the  most  important  collisional  process. 

From  the  point  of  view  of  radiative  transfer,  the  procedure  of 
'  orbvn  et  al.  [1970]  of  treating  the  radiative  transfer  problem  for 
a  single  molecular  vibrational-rotational  band  with  a  Doppler  line- 
shape  as  if  it  were  a  single  line  was  adopted.  As  they  point  out, 
this  procedure  greatly  simplifies  things  for  NLTE  conditions. 

Excited  N? ,  07 ,  and  0  Calculations 

Any  energy  deposited  in  the  *D  state  of  atomic  oxygen  will  undergo 
a  transfer,  at  least  to  some  extent,  to  N2  and  02 .  It  is  therefore 
necessary  for  one  to  determine  the  population  of  0(]D)  before  con¬ 
sidering  the  population  of  N2(v-1)  and  02(v=l). 

Rather  than  attempt  to  solve  time-dependent  equations  for  the  ex¬ 
cited  states,  it  was  decided  to  assume  that  an  equixibrium  had  been 

reached  between  the  production  and  loss  terms.  That  is,  lather  than 
solve 


dt 


Production  terms-loss  terms 


the  time  dependence  would  be  removed  by  assuming  quasi  steady  state,  i.e. 
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d£o ( lp) 3  _  0 

dt  U 

In  so  doing,  the  resulting  equation  to  be  solved  for  the  population  of 
O(lD)  is 

[0(]D)]  =  (kj [e][0( 3P) ]  +  k2[0t][e])  /  (qj[N2]  +  q2[02])  (52) 

The  two  production  terms  (k's)  are  obtained  from  the  program  EXCHNG  as 
described  in  a  previous  section.  The  loss  rate  terms  (q's)  which  are 
used  are  qj  =  9.0  x  10-11 
q2  =  4.0  x  10-1 1 

which  are  in  general  agreement  with  Young  et  al.  [1968]. 

Once  the  population  for  O^D)  has  been  determined,  it  is  possible 
to  compute  with  reasonable  accuracy  the  population  of  h2(v  =  1). 

In  making  this  calculation,  it  is  assumed  that  all  levels  with  v  >  2 
rapidly  decay  to  the  \  =  1  level.  Therefore,  one  need  only  be  concerned 
with  two  vibrational  levels;  v  =  0  and  v  =  1.  The  assumption  of  a  quasi 
steady  state  condition  leads  one  to  an  equation  for  N2(v=l),  i.e.  N^,  as 
follows : 

Of]  =  0cTV[M]  +  ke  +  k0[0(lD)]  +  k^^Cof]  +  k02[0f]  +  kE)[N2]/(qe  + 

qC02[C°2^  +  qN0[N°]  +  qN20[N2°]  +  qN0+i:N0+]  +  qH20[H20^  + 

+  qVT^M^  =  cJ0t0(3p)])  (53) 

Most  of  the  rate  constants  above  are  self  explanatory,  however,  a  few 
may  need  some  explanation. 

k^,y  -  production  rate  due  to  T— V  energy  exchange 
qVT  =  ^oss  ra^e  due  bo  V-T  energy  exchange 
kg  =  production  rate  due  to  secondary  electrons 
9e  =  loss  rate  due  to  ambient  electrons 
^  =  production  rate  attributable  to  electric  fields 
[m]  =  total  number  density  -  all  species 
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Figures  27  through  30  show  the  production  rate  kg  and  the  loss  rate 
qe  as  a  function  of  height.  The  first  two  figures  are  for  the  case  of 
an  electron  drizzle  but  no  aurora,  while  the  second  two  are  for  the 
case  of  an  IBC-III  aurora  five  minutes  after  commencement. 

The  same  assumptions  which  were  made  for  N2  were  also  made  in 
order  to  facilitate  computation  of  02(v=l),  i.e.  of.  The  quasi  steady 
state  condition  leads  to  an  equation  for  of  as  follows: 

[of]  "  <kTV[M]  +  ^  [4]  +  +  k0[0(lD)]  +  kE)[02]/(qe  +  + 

qo3[03l  +  qN  [n2]  +  qVT[M]  +  qAEX[o(3P)^)  (54) 

Here  the  production  and  loss  terms  take  on  the  same  meanings  as  out¬ 
lined  above.  The  values  of  the  constants,  needless  to  say,  are  no 
longer  the  same.  The  one  addition  is  the  q  term  which  accounts  for 

JL 

me  a  bum  uiLciuiduge  icatuuu  u  -r  u2  -►  0  +  02  as  discussed  by  Deggs 

[1971]. 

Radiative  Transfer  Considerations 

Contributions  to  the  broadening  of  an  absorption  line  are  in 
general  due  to  five  processes.  Each  of  these  processes  acts  to 
broaden  che  absorption  line.  Under  certain  conditions,  all  five  may 
be  important.  On  the  other  hand,  it  may  well  be  that  one  particular 
process  dominates  all  the  others.  The  five  processes  are  listed  be¬ 
low. 

1.  Doppler  Eroadening  -  the  result  of  molecular  or  atomic  motion 
(velocity  broadening) 

2.  Natural  Broadening  -  the  result  of  excited  states  having  a 
finite  lifetime. 

3.  Lorentz  Broadening  -  the  result  of  collisions  with  foreign 
gas  molecules . 

4.  Holtzmark  Broadening  -  the  result  of  collisions  with  absorbing 
molecules  of  the  same  kind. 


5.  Stark  Broadening  -  the  result  of  collisions  with  electrons  and 
ions . 
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Lorentz  ind  Holtsmark  Broadening  are  often  referred  to  as  pressure 
broadening.  Broadening  as  the  result  of  collisions  becomes  unimportant 
with  respect  to  Doppler  Broadening  at  altitudes  above  approximately  70  km 
[Kuhn  and  Lo*  1969].  Of  the  two  remaining  processes,  unless  the  op¬ 
tical  thickness  is  quite  large  (kQl  >  1000)  the  central  region  of  the 
line  will  dominate  the  absorption  process.  That  is  to  say,  the  edges 
of  the  absorption  line  may  be  neglected  and  the  entire  line  may  be  well 
approximated  with  a  purely  Doppler  lineshape.  If,  however,  the  absorbing 
gas  is  optically  very  thick,  then  only  in  the  edges  would  a  measurable 
amount  of  light  be  transmitted.  For  this  optically  thick  case.  Natural 
Broadening  must  be  considered. 

A  study  of  the  vibrational  rotational  bands  of  the  minor  constitu¬ 
ents  under  consideration  indicates,  however,  that  the  only  band  of  sig¬ 
nificant  optical  thickness  is  the  4.3  u  band  of  C02 .  Even  this  band 
does  not  approach  the  requirements  [ Mitchell  and  Zemansh  1971]  to  uti¬ 
lize  Natural  Broadening  as  a  lineshape  factor.  Therefore,  the  computa¬ 
tional  model  assumes  a  Doppler  lineshape  factor. 

The  spectral  absorption  coefficient  of  a  Doppler  broadened  line  is 
[ Mitchell  and  Zemcmsky  1971] 


r  2(v  -  v  )  _ ; 

k  =  k  exp  |  - - - - —  Jtn2 

v  o  L  Avd  J 
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where 


Av. 


2/2R£n2  pT 

C  "o  /M 


is  the  Doppler  breadth.  T  and  M  are  the  absolute  temperature  and  the 
molecular  weight  respectively.  kQ  is  the  maximum  absorption  coefficient 
which  occurs  at  v  =  vQ.  k0  may  be  found  from 


'\dv  •  ‘V®  kc4vD 

which  yields  the  expression 

k  =  _2_  /Zw 2  Xn2  82  N 
o  AvD  /  tt  8ir  gj  t 


(56) 


or 


(57) 


Here  the  parameter  g2/g}  is  the  ratio  of  the  statistical  weights  of 
the  upper  to  the  lower  state.  -  is  the  lifetime  of  molecule  in  the 
upper  state  and,  hence,  the  Einstein  A  coefficient  is  A  =  1/t.  The 
parameter  f  is  the  absorption  oscillator  strength.  Under  conditions 
when  hv^u  >>  kT,  f  and  the  band  strength  are  related  by  a  constant  as 
f Penner  1959] 


S,  =  2.3795  x  107  f 
lu 


(58) 


Ever,  though  the  optical  thickness  of  the  various  vibrational-rota¬ 
tional  bands  is  small  enough  to  adequately  describe  the  lineshape  by 
Doppler  broadening,  it  is  not  small  enough  to  neglect.  The  effect  of 
jptical  thickness  is  a  reduction  of  the  radiant  flux  being  transmitted 
through  the  medium.  Let  us  consider  the  radiative  exchange  process  in 
more  detail. 

For  a  parallel  beam  of  incident  intensity  Iy,  the  amount  of  energy 
absorbed  by  a  gaseous  column  of  length  Z  per  unit  time  is  [ Mitchell  and 
Zemansky  1961] 


F 


co  -k  Z 

f  I  (1  -  e  )  dv 
o  v 


(59) 


If  the  gas  is  Doppler  broadened,  this  expression  on  a  per  unit  length 
basis  becomes 


F 


OO 

/  I  (1  -  exp[-k  exp{- 
o  v  '  r  o  r  1 


2(  >  - 


v  ) 
o 


Av„ 


2 


JE3 


>])dv 
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If  we  define  w  = 


then  dw  = 


2(v  -  vo) 


2/Zn2 

Av„ 


dv 


Jt\ 12 


Now  the  line  will  absorb  only  near  its  center  frequency  vq  hence  we  may 
write 
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Av  -  2 

F  =  I^v0)  ZfZ  /  [l  ~  exp(-k  e-”  )]dw 
2  v£n2  0 


(61) 


Av. 


F-  Wr^=^0n~  +  -2r- - ]  (62, 

2  /In  2  2l/2  3!/3 


0  k  k  2. 

F  =  I  (v  Nf  [l  -  — +  — — 

o  °  me 


2.'/2  3 1  /3 

k  k  2 

Let  us  define  S(k  )  =  Tl - —  +  — 2 -  1 

O  .  T~~  r—  *  *  *  *  -* 

21/2  3I/J 


•  J 


(63) 


(64) 


A  factor  L(k^)  may  be  defined  as 


L(ko)  5  kQ  S(ko)  =  -  exp(-koe-w2)]dw  (65) 

/tt 

If  one  now  takes  the  derivative  of  L(kp)  with  respect  to  k  ,  it  will 

give  a  transmission  factor  for  an  optical  depth  kQ  as  [ Ivanov  and  Shcher¬ 
bakov  1965] 


M 


i(y 


dL(k0}  1  - 
~dk  =  7:  /_expC^2 

o  /m 


koe-»2]„ 


1W 


(66) 


In  a  similar  manner  one  may  obtain  transmission  factors  for  cases 
other  than  for  a  parallel  beam  continuum  source.  Specifically,  Nj  (k  ) 
and  N?(kQ)  are  used  in  this  report  and  represent  transmission  factors 
for  a  continuum  surface  source  and  a  doppler  lineshape  layer  source. 
Expressions  for  Nj (k  )  and  N2 (k  )  are  [ Deags  1972] 


2m  tt  /  2  k 

Nl(kQ)  =  /  d(f >/  cos0sin6Mi( — ^r)d0 

n  ^  COSQ 


(67) 


2m  m/2 


N2<ko>  '  /„  d*/0  «2fefe>  sl"9de 


(68) 


where 


M2<k,)  -  -p  IZ  Wl-*'2  -  k  e'”2]d» 

/tt 


(69) 


The  computationa  irocedure  used  to  determine  radiative  absorption 
was  adopted  from  'orbtr  et  al.  [1970].  In  essence,  the  effect  of  optical 
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thickness  is  accounted  for  by  multiplying  the  radiant  flux  (whether  it 
be  solar,  terrestrial,  or  computed  flux  radiated  from  some  small  atmo¬ 
spheric  volume)  by  one  of  the  three  transmission  factors.  These  fac¬ 
tors,  M^(x),  NjCt),  and  N^t),  are  dimensionless  and  are  computed  for 
all  vibrational-rotational  bands  (the  superscript  B  designates  the 
factor  for  a  band  rather  than  a  line)  by  means  of  a  best  fit  polynomial 
expression.  The  choice  of  factors  is  determined  by  the  kind  of  radi¬ 
ation  source.  That  is  to  say,  M^(t)  is  ’  ^sd  for  both  solar  flux  and 
for  the  particular  infrared  radiance  of  a  vibration-rotational  band 
which  we  are  trying  to  find.  Specifically,  M?(t)  is  used  when  the 
radiant  flux  is  in  the  form  of  a  parallel  beam.  On  the  other  hand, 
the  terrestrial  flux  is  essentially  a  black  body  source  which  is  2tt 
steradians  in  extent.  N^(t)  is  the  proper  factor  to  use  with  this 
type  of  source.  The  third  and  last  source  of  radiant  flux  is  that 
generally  categorized  as  other  parts  of  the  atmosphere.  Here  the 
source  of  radiation  has  a  Doppler  lineshape  profile  which  is  assumed 
to  be  exactly  that  of  the  absorbing  line.  In  addition,  this  source 
of  radiation  is  handled  in  the  computational  model  on  a  layer  by  layer 
basis.  The  proper  factor  for  this  type  of  source  is  ’^(i). 

As  previously  stated,  M?(t),  N?(t),  and  N^(x)  are  computed  by 
means  of  best  fit  polynomials  developed  by  Corbin  et  al.  [1970]. 

The  particular  polynomial  to  be  used  is  determined  by  (1)  the  optical 
thickness  x  and  (2)  whether  the  band  is  parallel  or  perpendicular 
transition.  Nevertheless,  all  bands  with  these  same  characteristics 
are  treated  in  the  same  manner.  At  first  this  may  be  viewed  with  some 
skepticism,  however,  in  light  of  the  optical  thinness  of  almost  all 
lines,  this  factor  turns  out  to  be  either  1  or  very  close  to  it. 

Therefore,  it  is  felt  that  little  error  is  introduced  by  this  method 
of  treating  the  problem. 

One^f urther  consideration  is  necessary  with  regard  to  the  factors 
Ml(T)»  N:(t),  and  N^t).  The  superscript  B  indicates  that  the  factor 
is  for  a  vibrational-rotational  band,  wnereas  in  general  it  is  written 
on  a  line  by  line  basis.  The  physics  of  transforming  these  factors  to 
a  band  basis  will  be  explained  below.  However,  some  of  the  technique 
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is  not  obvious  in  the  computer  program  because  it  is  included  in  the 
best  fit  polynomials  for  evaluating  M?(x),  Ni(x),  and  n|(t) . 

Since  the  relative  line-strength  of  each  line  in  a  viurational- 
rotational  band  varies  with  respect  to  the  other  lines,  so  will  the 
optical  thickness  of  each  line.  In  order  to  greatly  simplify  computa¬ 
tion,  only  one  value  of  optical  thickness  (Tq)  is  used  for  the  whole 
band.  This  mean;  that  only  one  computation  of  either  M?(x),  N?(t) ,  or 
N2(t)  is  necessary  and  that  that  value  will  then  be  on  a  band  basis. 

Thus  a  reasonable  methol  of  choosing  tq  is  all  that  is  outwardly  re¬ 
quired.  Inwardly,  on  the  other  hand,  the  polynomial  fit  method  shields 
from  view  the  process  used  to  obtain  it.  That  is  in  order  to  obtain 
M](t)  for  example,  the  product  tMi(t)  is  computed  for  each  line  and  then 
summed  for  all  lines  where  the  product  is  more  than  10“14  that  of  the 
strongest  line.  This  sum  is  divided  by  the  sum  of  the  optical  thickness 
t  of  all  these  lines  in  order  to  obtain  M?(t).  Mathematically  then, 


(70) 


where  i,  j,  and  k  are  over  lines  in  the  P,  Q,  and  R  branches. 

The  optical  thickness  tq  is  assigned  to  be  that  of  the  strongest 
line.  This  value  is  determined  in  the  following  manner.  The  line  in¬ 
tensity  is  proportional  to  the  fractional  part  of  molecules  in  the  ini¬ 
tial  state  [ Goody  1964]. 

n(J")  gJ"  .  Er(J"), 

n  Qr(t)'  exp  L-  kT  ]  (71) 

where  Qr(T)  is  the  rotational  partition  function  at  the  temperature  T 
which  is  given  by 

"  E  (J) 

Qr(T)  =  j=o  8j  exp  (72) 

The  term  gj  is  the  multiplicity  of  the  quantum  state  J.  E(J)  is  the 
rotational  energy  level  of  the  molecule.  For  example,  the^otacional 
energy  levels  of  a  rigid  linear  molecule  are  given  by  [Goody  1964] 

Er(J)  =  Bv  he  J(J  +  1)  (73) 
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The  quantium  number  J"  for  the  strongest  line  may  be  found  by  the  max¬ 
imization  process  of  setting  the  derivative  of  n(J)/n  with  respect  to 

J  equal  to  zero  and  solving.  The  result  is  that  the  maximum  line  in¬ 
tensity  occurs  at 

_  /  kT  , 

'  /  2B  he  "  **  (74) 

v 


Thus,  by  using  this  value  of  J"  one  may  determine  the  relative 
strength  of  the  strongest  line  by  computing  n(J")/n.  The  value  t 
may  then  be  found  using  equation  (57) 


e^/mc  /  tt"1  /m 

2.3795  x  10 7  /  2Rj/  T  Nl'VS’SM  (75) 

where  S-SM  is  the  band  strength  (S^)  mUi  tiplied  by  the  relative 
strength  of  the  strongest  line  (n(J")/n). 


NLTE  Vibrational  Populations 

In  an  effort  to  compute  the  vibrational  populations  of  the  various 

species,  let  us  first  list  the  possible  excitation  or  de-excitation 
mechanisms . 

Excitation  (production  terms) : 

a)  collisional  T-V  excitation 

b)  collisional  V-V  excitation 

c)  radiative  absorption  of  solar  flux 

d)  radiative  absorption  of  terrestrial  flux 

e)  radiative  absorption  of  flux  from  other  levels  of  the 
atmosphere 

f)  chemiluminescence 

g)  spontaneous  radiative  transition  from  a  higher  state 
De-Excitation  (loss  terms) : 

a)  collisional  V-T  de-excitation 

b)  collisional  V-V  de-excitation 

c)  spontaneous  radiative  de— excitation 

With  these  processes  in  mind  it  is  possible  to  write  a  general 
differential  equation  for  the  time  rate  of  change  of  the  number  density 
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of  the  upper  state  dNu/dt  in  terms  of  them 

dNu/dt  =  production  terms  -  loss  terms 

£  4= 

As  in  the  computation  for  N2  and  O2  we  will  assume  a  case  of  quasi 
steady  state.  That  is,  we  will  assume  that  dN^/dt  =  0  and  then  solve 
the  resulting  equation  for  N^,  the  vibrational  population  of  the  state 
under  consideration. 

Since  all  of  the  terms  are  not  present  in  every  case,  let  us  con¬ 
sider  a  typical  example  of  the  computation  of  a  vibrational  population. 

Suppose  we  consider  the  population  and  depopulation  of  the  tt  level  of 

u 

the  V2  bending  mode  for  CO2  (see  Figure  32)  .  Here  one  would  write : 
C02(iru)  =  ( (kTV  +  GUp*N?(T)  +  2irG-Neff*A1  +  0^(1))  ([C02])  + 

[C02ag)](A2  +  q^)  +  [C02(Ag)](A3  +  q^)  + 

[C°2(Eg(v1))](A7  +  qVT)}/(qVT  +  A3)  (76) 

In  this  expression  and  q^  are  the  proper  T-V  and  V-T  collisional 
energy  exchange  rate  coefficients. 

G^  =  G*terrestrial  flux*2ir  (sec-1) 

GDOWN  =  G*solar  flux  (sec-1) 

2 

where  G  =  f  A  2  =  3.72  x  10-34  SA  2  (cm2«p) 

mc^  o  o 

GUP  anc^  GDOWN  are  derived  from  the  general  expression  of  the  rate  at 
which  a  single  molecule  absorps  and  radiates  infrared  radiation.  Equa¬ 
tion  (59)  gives  their  rate  as 

G1  =  Q  Iv  (1  -  e'kv)  dv  (77) 

on  a  per  unit  length  basis.  To  transform  this  expression  to  a  wave¬ 
length  basis,  we  must  recall  that  dA  =  A2/c  dv.  Then  one  can  "rite 
this  rate  as 


G1  =  f°  1(A)  —  (1  -  e"k*)dA 
'  o  c 


(78) 
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which  upon  integration  (since  N  =  1)  yields 


'1  . 


1(V  4  f]  -  «»„>  „c 


TTe 


fX  2 
o 


(79) 


provided  that  I(Aq)  has  not  diminished.  That  is  to  say,  we  will  account 
for  absorption  by  multiplying  by  the  appropriate  transmission  factor  as 
in  equation  (76)  rather  than  concerning  ourselves  with  it  at  this  point. 
Thus  it  can  be  seen  that  and  are  just  special  cases  of  G> . 

Neff  is  the  effective  number  of  molecules  contributing  to  the 
radiative  flux  from  other  parts  of  the  atmosphere,  in  a  given  vibra¬ 


tional-rotational  band,  multiplied  by  the  factor  N](t).  N  has  the 
units  photons  cm"'  p"1 .  AT  is  the  Einstein  A  coefficient  for  spon¬ 
taneous  -mission.  Note  that  all  of  the  factors  up  to  this  point  are 
multiplied  by  [C02],  which  gives  the  resulting  units  cm-3  sec-1  as  de¬ 
sired.  The  term  [C02 (1+) ]*(A2  +  q^)  and  the  two  succeeding  terms 
give  the  production  rate  of  the  desired  level,  C02  (t^)  ,  (see  Figure  32 
for  an  energy  level  diagram)  from  higher  levels.  This  production  may 
be  accomplished  via  spontaneous  emissions  (A2)  or  through  collisional 
energy  loss  (q^)  .  The  loss  terms  for  the  ttu  level  of  the  v2  bending 
mode  of  CO 2  are  obvious  in  light  of  the  previous  discussion. 

It  should  be  pointed  out  that  all  of  the  equations  such  as  this 
one  are  solved  by  means  of  an  iterative  process.  The  initial  values 
for  the  number  density  of  the  various  energy  levels  used  in  this  pro¬ 
cess  are  those  as  determined  from  a  Boltzman  distribution. 


Chemiluminescence 

Infrared  chemiluminescence  may  be  defined  as  infrared  energy  radi 
ated  as  the  result  of  a  chemical  reaction  having  taken  place.  This 
statement  implies  that  at  least  one  of  the  product  molecules  of  the 
chemical  reaction  is  left  in  an  excited  vibrational-rotational  state 
which  subseuqently  radiates  in  the  infrared.  This  means,  for  example, 
thac  when  the  following  chemical  reaction  occurs 

+  +  ± 

0  +  N2  -*•  (NO  )  +  N 


(80) 


78 


Figure  33. 
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the  product  molecules  have  less  internal  energy  in  their  respective 
ground  states  than  do  the  reacting  molecules .  In  order  to  satisfy  the 
law  of  conservation  of  energy,  part  of  this  additional  energy  may  go 
into  vibrational-rotational  energy  and  the  remaining  part  into  the 
translational  energy  of  the  product  molecules.  In  the  case  of  the 
above  reaction,  we  have  assumed  that  one  quantum  of  vibrational-rota¬ 
tional  energy  is  obtained  per  reaction.  Further,  we  assume  that  1/3  of 
this  goes  into  the  v=l  vibrational  level  and  1/3  into  the  v=2  level. 
Since  the  v=2  level  is  approximately  twice  the  energy  of  the  v=l  level, 
the  total  of  the  two  is  1  quantum  of  energy. 

The  production  rate  of  NO  from  this  reaction  is  obtained  f^om  the 
program  EXCHNG.  This  production  rate  multiplied  by  the  proper  factor 
(1/3  for  example)  is  then  included  in  as  an  additional  production  term 
in  the  equation  used  to  compute  vibiational  populations. 

Two  other  chemiluminescent  reactions  are  considered.  They  are 


N  +  02  NO*  +  0 
H  +  03  ->  OH  +  02 


(81) 

(82) 


In  the  former,  atomic  nitrogen  may  be  considered  to  be  in  either  the 
N(t*S)  or  the  N(^D)  states.  We  have  taken  the  production  rates  given  by 
Ghosh  [1968]  and  modified  them  by  multiplying  by  0.02.  This  is  done 
for  both  the  v=l  and  the  v-2  vibrational  levels.  This  essentially  limits 
us  to  considering  only  the  N(l+S)  reaction.  Less  is  known  regarding  the 
N(2D)  reaction  and  essentially  no  accounting  of  it  is  made  in  this  pro¬ 
gram. 

The  hydroxyl  radical  OH  is  different  than  all  other  species  in  the 
program  in  that  the  only  important  nighttime  excitation  mechanism  is 
through  chemiluminescence.  This  is  the  result  of  the  nonexistence  of 
a  near  resonant  V-V  exchange  mechanism  and  of  the  fact  that  the  ter¬ 
restrial  flux  in  this  spectral  region  is  relatively  very  small.  In  order 
to  compute  the  vibrational  populations  of  the  upper  vibrational  levels 
of  OH,  it  was  assumed  that  the  above  chemical  reaction  was  the  only  one 
of  any  consequence.  Further,  it  should  be  pointed  out  that  this  reac¬ 
tion  is  capable  of  populating  only  the  levels  up  to  v=9 .  A  production 
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rate  of  5.69  x  10  ^  cm^  molecule  '  sec  ^  was  used  for  the  v=9  vibra¬ 
tional  level  and  a  smoothed  version  of  the  relative  formation  rates  of 
Charters  et  al.  ri9?i]  was  adopted  as  a  means  of  determining  lower  vi¬ 
brational  level  production  rates. 

PROGRAM  BCKGND  (IR  BAND  MODEL  PROGRAM)  RESULTS 

It  is  the  purpose  of  this  section  to  present,  in  graphical  form, 
some  of  the  results  of  the  program  BCKGND.  It  will  be  recalled  that 
this  program  has  the  capability  of  computing  the  predicted  radiance  in 
watts  cm  ^  ster  ^  on  a  band  basis.  It  also  has  the  versatility  of  com¬ 
puting  this  radiance  for  the  angles  0°,  20°,  40°,  60°,  70°,  80°,  and 
90°,  as  measured  from  the  zenith,  as  well  as  for  the  limb  view  case. 
Needless  to  say,  this  Is  a  large  amount  of  data.  In  an  effort  to  re¬ 
duce  the  total  amount  of  data  presented  In  this  report,  and  yet,  to 
keep  it  representative  and  meaningful,  it  was  decided  to  present  for 
the  most  part  only  data  In  which  the  zenith  angle  is  0°.  This  is  to 
say,  only  vertical  viewing  data  will  be  presented. 

The  predicted  radiance,  for  the  vertical  viewing  case,  is  equiva¬ 
lent  to  predicting  what  radiance  a  forward  looking  radiometer  on  a  rocket 
would  observe.  It  is  still  possible  to  see  the  effects  of  optical  thick¬ 
ness  for  this  ease  as  compared  to  an  optically  thin  medium.  It  is  also 
possible  to  determine  how  the  increase  in  electron  flux  affects  the  pre¬ 
dicted  radiance  and  in  particular  to  observe  which  bands  are  enhanced 
by  the  auroral  activity.  The  same  thing  may  be  said  for  the  case  in 
which  the  energy  input  is  due  to  an  electric  field. 

Let  us  first  consider  the  case  in  which  there  is  no  aurora,  no 
electric  fields,  but  in  which  an  electron  drizzle  (as  previously  de¬ 
scribed)  exists.  The  predicted  radiance  for  vertical  viewing  is  shown 
in  Figures  34  through  44  for  the  molecules  CH4,  C02,  H20  vapor,  NO, 

N20,  O3,  NO  ,  OH,  and  CO.  The  designation  NLTE=3  means  that  nonthermo¬ 
dynamic  equilibrium  conditions  exist.  In  this  case  the  3  means  night¬ 
time  conditions.  The  designation  aurora  =  0  means  no  auroral  conditions 
exist. 
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The  radiance  from  CO2  is  shown  in  two  graphs.  Figure  36  indicates 
that  the  plot  is  of  the  15  p  bands.  This  means  that  all  of  the  bands 
near  15  p  have  summed  together  and  are  being  plotted  as  one  collective 
band . 

The  two  plots  of  NO  shown  in  Figures  38  ard  39  are  identical  plots. 
The  region  from  60-150  km  is  redrawn  on  an  expanded  scale  to  facilitate 
comparison. 

Only  two  of  the  Av  =  2  bands  of  OH  are  shown  in  Figure  43.  This  is 
done  simply  as  a  means  of  reducing  the  total  number  of  bands  to  be 
plotted.  One  could,  if  desired,  include  all  Av  =  2  bands  as  well  as 
the  Av  =  1,  etc.  This  is  done  in  Figure  78. 

As  a  means  of  determining  the  effect  of  the  electron  drizzle  on 
the  radiance  of  the  various  bands,  let  us  compare  Figures  35  and  40  with 
Figures  45  and  46.  The  latter  two  figures  are  plots  of  the  predicted 
radiance  of  two  of  the  CO2  bands  and  the  N2O  bands  with  no  aurora,  no 
electric  fields,  and  no  electron  drizzle.  It  can  be  seen  that  during 
the  election  drizzle  both  CC2  bands  are  somewhat  enhanced  but  that  there 
is  only  a  very  slight  enhancement  of  the  4.5  p  band  of  N20.  Neither  the 
7.78  p  band  nor  the  17.0  p  band  of  N2O  is  enhanced.  By  referring  to  the 
energy  level  diagrams  of  Figures  32  and  33,  one  can  see  that  the 
cause  of  the  selective  enhancement  is  the  near  resonant  V-V  energy  ex¬ 
change  with  N2 .  That  is  to  say,  the  electron  drizzle  tends  to  increase 
the  vibrational  energy  of  N2  which  in  turn  exchanges  part  of  this  energy 
with  selected  levels  of  CO2  and  N2O.  These  levels  in  turn  radiate,  and 
having  an  enhanced  population,  give  rise  to  an  increased  radiance. 

Next  let  us  see  what  effect  an  IBC-II  aurora  has  on  the  predicted 
radiance.  Figures  47  through  55  are  plots  of  the  predicted  radiance 
for  the  case  of  an  IBC-II  aurora  five  minutes  after  commencement  plus 
the  background  electron  drizzle  but  with  no  electric  field.  These  plots 
represent  the  radiance  viewing  vertically  upward.  The  designation  AURO¬ 
RA  =  1  simply  means  that  an  IBC-II  aurora  exists.  The  auroral  region 
was  taken  to  be  from  80  km  of  -.ltitude  to  150  km.  The  integration  limits 
in  determining  the  radiance  were  taken  to  be  these  value  .  Theiefore, 
it  should  be  pointed  out  that  this  represents  the  total  radiance  (aurora 
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Figure  53.  Height  profile  of  the  band  radiance  of  NO  five  minutes  after 


profile  of  the  band  radiance  of  two  of  the  bands  of  OH  five  minutes  after 
an  IBC-II  aurora. 
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+  background)  for  every  case  except  NO.  The  reason  for  NO  being  an  ex¬ 
ception  is  that  it  is  the  only  infrared  active  gas  with  significant 
number  densities  at  altitudes  greater  than  150  km.  Hence,  in  the  case 
of  NO  the  radiance  plotted  in  Figure  50  represents  only  the  auroral 
contribution;  i.e.,  the  background  contribution  from  molecules  at  al¬ 
titudes  greater  than  150  km  is  not  included. 

It  can  be  seen  that  the  aurorally  enhanced  bands  include  C02  at 
4.3  p,  10.4  u  (and  some  of  the  15  p  band  group),  NO  at  5.32  p,  N20  at 
4.5  p,  and  all  of  the  N0+  bands.  The  molecule  CO  has  a  near  resonant 
V-V  exchange  with  N2  but  the  number  densities  of  this  gas  at  auroral 
altitudes  is  not  sufficient  to  cause  significant  auroral  enhancement. 

Of  the  remaining  gases,  CH4,  H20,  and  03  have  a  near  resonant  V-V  ex¬ 
change  with  02  rather  than  N2 ,  and  OH  has  no  near  resonant  V-V  exchange 
with  either  of  these  molecules.  The  much  smaller  cross  sections  for 
an  electron-02  collision  (as  compared  with  N2)  make  the  auroral  en¬ 
hancement  of  the  latter  gases  insignificant. 

Next  let  us  increase  the  auroral  activity  to  that  of  an  IBC-III 
aurora  (AURORA  =  3)  and  compare  the  predicted  radiance  with  previous 
results.  Figures  56  through  65  show  the  predicted  radiance  of  an  IBC- 
III  aurora  five  minutes  after  commencement  with  the  background  electron 
drizzle  but  r;ith  no  electric  fields.  As  before,  the  NLTE  *  3  designa¬ 
tion  means  nighttime  NLTE  conditions.  Again,  one  may  note  that  there 
is  a  further  enhancement  of  C02  bands,  the  NO  band  at  5.32  p,  the  N20 
band  at  4.5  p,  and  all  of  the  N0+  bands. 

For  the  next  case  let  us  go  back  to  quiescent  conditions  but  add 
an  electric  field.  Figures  66  through  74  show  the  predicted  radiance 
when  no  aurora  or  electron  drizzle  is  occurring  but  where  an  electric 
field  of  0.1  volts/meter  is  present.  The  electric  field,  as  described 
in  a  previous  section,  leads  to  a  vibrational  excitation  of  both  N2 
and  02.  Again  the  near  resonant  V-V  exchange  process  becomes  important 
and  leads  to  an  enhancement  of  various  bands.  In  the  case  of  an  elec¬ 
tric  field  as  the  disturbing  element,  the  02(v=l)  level  is  more  heavily 
populated  leading  to  some  enhancement  of  CH4  at  6.55  p,  H20,  and  03  in 
addition  to  the  enhancement  of  C02,  NO,  N20,  and  N0+  as  before.  One 
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Figure  65.  Height  profile  of  the  band  radiance  of  CO  five  minutes  after  commencement  of  an 
IBC-Itl  aurora. 
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Figure  72.  Height  profile  of  the  band  radiance  of  NO  for  the  case  when  the  only  source  of 
ionospheric  disturbance  is  an  electric  field  of  .0.1  v/m. 
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might  note  that  the  altitude  regime  for  enhancement  is  somewhat  differ¬ 
ent  than  that  of  the  aurorally  disturbed  cases.  This  is  because  the 
electric  field  excitation  process  cannot  produce  02(v=l)  or  N2 (v=2)  at 
altitudes  where  tne  mean  free  path  is  below  some  threshold  length.  That 
is  to  say,  the  electric  field  excitation  process  has  a  threshold  alti¬ 
tude  of  somewhere  near  100-113  km  or  so. 

As  a  comparison,  Figures  75  and  76  show  the  predicted  radiance  of 
C02  at  4.3  u  and  10.4  y  and  the  bands  of  NO+  when  the  electric  field  has 
been  increased  to  0.2  volts/meter. 

Much  has  been  said  abo  ■  the  near  resonant  V-V  exchange  process 
being  the  means  of  enhancing  certain  vibration  populations  of  the  infra¬ 
red  active  specxes.  The  amount  of  energy  exchange  to  the  minor  constit¬ 
uents,  therefore,  is  proportional  t(  the  vibrational  population  of  the 
v-1  state  of  N2  (or  02  as  the  case  may  be)  .  Figure  77  is  a  plot  of  the 
ratio  of  the  predicted  v=l  level  nopulation  of  N2  to  the  Boltzman  popu¬ 
lation  of  that  level.  Five  different  cases  are  presented  in  this  plot. 
They  are  (1)  totally  quiescent  conditions,  (2)  that  of  only  a  background 
electron  drizzle,  (3)  IBC-II  aurora  five  minutes  after  commencement  with 
electron  drizzle,  (4)  IBC-III  aurora  five  minutes  after  commencement 
with  electron  drizzle,  and  (5)  that  with  an  electric  field  of  0.1  volts/ 
meter  but  no  aurora  or  electron  drizzle.  Needless  to  say,  the  N  (v=l) 
population  significantly  increases  with  auroral  activity.  One  might 
also  note  the  altitude  dependence  of  this  population  with  the  excita¬ 
tion  mechanism. 

SPECTRUM  COMPUTATION  PROGRAM  (PROGRAM  SPCTRA)  PHYSICS 

The  data  presented  in  the  previous  section  was  on  a  band-by-band 
basis.  This  is  informative  if  one  is  interested  in  only  one  particular 
band.  Very  often,  however,  it  is  of  interest  to  know  the  predicted  in¬ 
frared  radiance  throughout  some  spectral  region.  In  chis  case  one  would 
need  to  know  the  radiance  contribution  of  each  and  every  band  which  radi¬ 
ates  in  that  particular  region. 

It  is  the  purpose  of  program  SPCTRA  to  commute  a  predicted  radiance 
throughout  the  spectral  region  1-25  L  .  This  is  lone  by  using  the  output 
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of  program  BCKGND  as  an  input  to  this  program.  That  is,  each  of  the 
predicted  band  radiances  for  each  of  the  minor  constituents  is  redis¬ 
tributed  on  a  line  by  line  basis  All  of  the  myriad  of  lines  from  all 
of  the  bands  of  all  of  the  minor  constituents  are  then  summed  together 
and  printed  out  every  0.1  u  from  1.0  u  to  25.0  p.  Three  different  view- 
ing  teometries  are  considered.  The  observer  may  be  looking  vertically 
upwards,  horizontally,  or  he  may  be  observing  in  limb  view  fashion. 


Line  Position  Determination 

The  rotational  energy  levels  for  a  rigid  linear  molecule  are  given 
as 


E  =  he  F( J) 

(83) 

where 

F(J)  =  B(J  +  1)J 

(84) 

B  is  the  rotational  constant  ft  that  particular  vibrational  level  under 
consideration  and  J  is  the  lotational  quantum  number.  For  this  type  of 

molecule  (if  we  ignore  any  selection  rules  associated  with  symmetry 

properties)  the  line  positions,  in  wave  numbers,  for  a  vibrational-rota¬ 
tional  band  are  given  as  [Gocdy  1964] 

P(J")  =  vq/c  -  (B '  -  B")J"  +  (B*  -  B") J"2  (85) 

Q(J")  =  vq/c  +  (B’  -  B") J"  +  (B’  -  B") J"2  (86) 

R(J")  =  vq/c  +  2B '  +  (3B '  -  B") J"  +  (B'  -  B")J"2  (87) 

where  B'  and  I'  (B"  and  J")  are  the  rotational  constant  and  the  rota¬ 
tional  quantum  number  in  the  upper  (lower)  energy  state  respectively 
apP  Vo  t^ie  ^re9Ut  ncy  of  the  band  center.  The  letters  P,  Q,  and  R 
de.xnt  the  branch  of  the  band.  That  is,  if  AJ  =  J'  -  J"  is  taken  as 

AJ  -  -1  P-Branch 

AJ  =  0  Q-Branch 

AJ  =  +1  R-Branch 

In  addition  to  the  above,  the  vibrational  angular  momentum  quantum 
number  plays  an  important  role.  Transitions  with  regard  to  this  quantum 
number  fall  into  one  of  two  categories:  a£  =  0  or  A£  =  ±1 . 
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If  A-C  =  0  the  transition  is  called  a  parallel  band.  If  in  addition 
^  =  0  in  both  the  upper  and  lower  states  it  is  called  a  l -l  transition 
and  has  r.c>  Q-branch.  Both  the  P-branch  and  the  R-branch  turn  out  to  be 
regular  arrays  of  lines  with  fairly  constant  spacing.  A  gap  occurs  at 
the  band  center  frequency  v  .  The  R-branch  lies  on  the  high  frequency 
side  of  and  the  P-branch  on  the  low  frequency  side.  The  intensities 
of  thest  two  branches  are  not  symmetrical  due  to  the  fact  that  the 
first  line  on  either  side  of  (R(0)  and  P(l))  have  different  initial 
states;  i.e.,  I  is  different,  and  hence  have  different  populations  and 
also  the  fact  that  the  intensity  varies  with  the  fourth  power  of  the 
frequency.  If  t  t  0  in  the  upper  and  lower  states,  a  parallel  band 
still  results  but  it  is  no  longer  a  l-l  transition.  In  this  ~ase  a 
weak  Q-branch  will  be  exhibited.  No  fundamentals  are  of  this  type 
[tizrzberg  1945]. 

If  =  ±1,  the  transition  is  called  a  perpendicular  band.  Bands 
of  this  type  exhibit  a  Q-branch  which  is  stronger  than  either  the  P- 
branch  or  the-  R-branch. 

Of  the  molecules  examined  in  this  report,  all  are  linear  with  the 
exception  of  CH4 ,  HpO,  and  O3.  Therefore,  the  line  positions  of  all 
but  these  may  be  computed  using  the  method  outlined  earlier  in  this 
section.  The  remaining  are  dealt  with  in  one  of  two  ways.  Either  the 
line  positions  are  obtained  from  experimental  measurement  and  used  as 
an  input  parameter,  or  individual  line  positions  are  not  computed.  In 
the  latter  case,  a  technique  of  computing  the  total  radiance  of  all 
lines  in  certain  spectral  intervals  is  used. 


Line  Strength  Determination 

In  a  previous  section  it  was  stated  that  the  line  intensity  is 
proportional  to  the  fractional  part  of  molecules  in  the  initial  state 
which  is  given  by 


n(J") 


n 


8J" 

Q  (T) 


t  U  ) 

exp 


(88) 


Therefore  the  individual  line  intensities  are  proportional  to 

he  B 

exp  [- 


kT 


J(J  +  D] 


(89) 
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multiplied  by  a  degeneracy  factor.  By  computing  this  intensity  i  -r 
each  line  and  summing  over  all  lines  a  total  relative  line  strength  for 
all  lines  of  all  branches  may  be  computed. 

The  spectral  absorption  coefficient  i  for  each  line  may  be  found 

I-j 

in  the  same  manner  as  equation  (75)  with  the  exception  that  a  line 
strength  and  position  must  be  used  rather  than  a  band  strength  and  posi¬ 
tion.  That  is, 

t  =  4.8315  x  10-18  /mTt  N  St  (90) 

l*  L  L 

where  for  example  if  a  line  in  the  P-branch  is  under  consideration 
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Here  PS  =  relative  line  strength 

'PS  =  total  relative  line  strength 
i 

Sg  =  band  strength 

The  predicted  radiance  for  a  particular  line  is  computed  as  fol¬ 
lows.  First,  all  of  the  relative  line  strengths  are  corrected  for  the 
effective  optical  thickness.  The  corrected  line  strengths  are  then 
summed  over  all  lines  and  all  branches  of  the  band.  The  predicted  line 
radiance  R,  is  then  found  as 

TLS  R 
i  1 


where  LS  =  the  corrected  relative  line  strength 

ELS  =  the  total  corrected  relative  line  strength 
i 

R  =  the  band  radiance  (watts-cm-2  -ster-1)  from  the  program 
BCKGND 


In  the  case  of  nonlinear  molecules  either  the  individual  line 
strengths  or  representative  line  strengths  throughout  the  band  inter¬ 
vals  are  given  as  input  variables.  That  is  to  say,  actual  computation 
of  individual  line  strengths  is  not  done  except  in  the  case  of  linear 
molecules . 
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Spectrum  Computation 

Once  the  computation  of  all  of  the  individual  line  positions  and 
their  predicted  radiance  has  been  completed  for  all  of  the  bands,  it 
is  a  simple  matter  to  compute  the  total  spectrum.  The  spectral  region 
from  1  p  to  25  u  is  essentially  broken  up  into  intervals  of  0.1  p.  The 
radiance  contributions  for  each  of  these  intervals  is  summed  over  all 
lines  of  all  bands  of  all  molecules.  As  mentioned  earlier,  this  is 
done  for  three  different  viewing  geomei  'ies  . 

RESULTS  OF  PROGRAM  SPCTRA 

The  output  of  program  SPCTRA,  which  computes  the  predicted  radi¬ 
ance  spectrum  from  1.0  p  to  25.0  p,  takes  three  different  geometries 
into  consideration.  Tuese  geometries  are  the  vertical  viewing  case, 
the  horizontal  viewing  case,  and  the  limb  viewing  case.  Each  of 
these  cases  is  done  for  the  altitudes  from  60  km  to  120  km  in  5-km 
increments.  As  in  a  previous  section,  it  is  not  necessary  to  present 
all  of  this  data  in  order  to  gain  some  insight  into  the  radiance  changes 
due  to  auroral  activity.  Instead,  two  representative  cases  are  illus¬ 
trated.  To  be  consistent  with  the  data  presented  from  program  BCKGND, 
the  vertical  viewing  case  was  chosen;  and  since  the  auroral  activity 

was  not  modeled  below  80  km,  this  altitude  was  chosen  as  the  viewing 
altitude. 

Figure  78  shows  the  predicted  radiance  spectrum  for  an  instrument 
at  80  km,  viewing  in  the  vertical  direction,  under  conditions  of  an 
electron  drizzle  but  no  aurora  and  also  for  an  IBC-III  aurora  five 
minutes  after  commencement.  The  figure  illustrates  how  the  predicted 
radiance  nanges  when  auroral  conditions  prevail.  The  electric  field 
intensity  is  assumed  to  be  zero • 

One  may  note  that  the  enhancement  due  to  the  auroral  energy  input 
apperas  in  the  bands  of  C02  at  4.3  p  and  10.4  p,  the  NO  band  at  5.3  p, 
the  N20  band  at  4.5  p,  and  the  N0+  bands  at  4.32  p,  4.27  p,  and  2.15  p 
as  predicted  by  the  band  model  program.  In  addition  there  is  an  en¬ 
hancement  at  13.9  p.  -This  raay  explained  as  follows. 


X  (microns) 


Figure  78.  The  solid  line  shov/s  the  rac. iance  of  the  infrared  spectrum 
from  1-25  u  under  conditions  of  an  electron  drizzle  only.  The  dashed 
line  indicates  the  enhancement  due  to  an  IBC-III  aurora.  The  predicted 
radiance  for  tiie  auroral  case  is  for  f.ve  minutes  after  auroral  com¬ 
mencement.  Doth  cases  are  for  nighttime  conditions. 
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Figure  36  is  a  plot  of  the  15  p  bands  of  C02  .  This  is  a  collection 
of  seven  different  bands  all  grouped  around  15  u.  Of  these,  the  most 
aurorally  enhanced  is  the  band  at  13.9  p.  The  total  radiance  of  this 
band,  however,  is  considerably  less  than  one  of  the  14.98  p  bands 
during  quiescent  conditions.  The  radiance  is  still  slightly  smaller 
than  this  14.98  p  band,  even  though  greatly  enhanced  by  the  IBC-III 
aurora.  One  would,  therefore,  expect  to  see  an  enhancement  at  13.9  p 
on  the  spectrum  plot  whereas  this  enhancement  is  hidden  from  view  in 
S^oiiping  all  of  the  bands  around  15.0  p  together. 

Plots  such  as  those  in  Figure  78  are  helpful  in  the  design  of 
scanning  radiometers  such  as  those  using  circularly  variable  filters. 
They  are  also  helpful  in  determining  filter  requirements  for  single 
channel  radiometers. 
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